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ABSTRACT 

This  thesis  presents  an  analytical  technique  for  the 
design  of  compensation  for  multi-loop  servomechanisms.  A 
completely  general  multi-loop  system  is  represented  in  a 
standard  block  diagram.   For  this  system  a  standard  deter- 
minant is  formed  which  uniquely  represents  each  block,  loop 
and  path  in  the  block  diagram.   Based  on  this  determinant, 
the  effect  of  all  possible  compensators  on  the  coefficients 
of  the  system  characteristic  equation  is  found.   This  infor- 
mation is  presented  in  a  Location  Matrix.  A  model  character- 
istic equation  is  synthesized  from  the  static  and  dynamic 
performance  specifications.   Comparison  of  the  uncompensated 
characteristic  equation  to  the  proper  model  equation  shows 
which  coefficients  must  be  altered.   The  Location  Matrix 
will  show  all  the  compensator-path  combinations  which  will 
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affect  the  desired  change.   The  method  requires  no  intuition, 
and  so  can  be  applied  by  persons  without  extensive  training 
in  the  servomechanism  field.  A  step-by-step  procedure  for 
application  of  this  technique  is  included  as  a  summary. 
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OBJECT 

The  object  of  this  thesis  is  to  present  an  analytical 
technique  for  the  design  of  multi-loop  servomechanism 
compensation.   This  technique  is  set  forth  under  the 
following  topics: 

(a)  Formation  of  the  Characteristic  Equation. 

(b)  The  Effect  of  Performance  Specifications 
on  the  Coefficients  of  the  Characteristic 
Equation. 

(c)  The  Effect  of  Compensation  on  the  Coef- 
ficients of  the  Characteristic  Equation. 

(d)  The  Choice  of  Compensation. 


Chapter  1 

INTRODUCTION 

and 

METHOD  OF  APPROACH 

lol  Introduction.   The  problem  of  compensation  design  for  a 
multi-loop  servomechanlsm  Is  not  a  new  one.   Indeed,  most 
standard  servomechanisms  or  control  systems  texts  give  It 
some  mention.   Most  such  works  treat  the  problem  as  an 
extension  of  the  slngle-lbop  problem,  and  attempt  to  apply 
single-loop  concepts  to  Its  solution. 

The  penalties  of  such  an  approach  are  numerous  and 
serious.   The  first  such  penalty  arises  from  the  technique 
of  getting  the  system  Into  a  single  loop  configuration.   This 
involves  a  block  diagram  reduction  or  some  similar  procedure. 
As  a  result  of  this  reduction,  the  Identities  of  the  original 
blocks,  and  the  physical  components  they  represent,  are  lost. 
Far  more  serious  than  this,  in  a  compensation  design  problem, 
is  the  fact  that  the  information  about  the  arrangement  of  the 
components  is  lost.   The  analysis  of  this  reduced  system  would 
show  algebraically  what  kind  of  compensation  was  needed.   It 
would  be  impossible  to  transform  this  information  into  a  given 
compensator  fed  around  certain  main  path  elements.   The  whole 
sense  of  element  and  path  configuration  no  longer  has  any 
correlation  with  the  algebra. 

There  are  two  ways  out  of  this  dilemma o   The  designer 
can  insert  all  possible  compensators,  in  all  usable  paths, 


and  re-do  the  analysis  until  he  finds  one  combination  that 
gives  the  desired  result.   Obviously  this  is  a  monumental 
task  for  any  system  of  reasonable  complexity.   There  are 
graphical  techniques  which  eliminate  some  of  this  work,  but 
they  still  involve  trial  and  error.   The  alternative  is  the 
intuitive  approach.   Here  the  designer  must  have  such  vast 
experience  that  he  generally  knows  the  effect  of  all  com- 
pensators, and  so  can  eliminate  from  consideration  many 
which  could  not  work. 

In  this  thesis  an  analytical  technique  is  presentee'  which 
requires  no  intuition.   It  examines  the  effect  of  all  possible 
compensators  without  the  accompanying  drudgery. 

1 . 2   Method  of  Approach 

A  general  multi-loop  control  system,  consisting  of  all 
possible  direct,  fredback  and  feed  forward  paths,  is  formulated. 
Its  components  are  presumed  to  be  describable  by  ordinary, 
linear  integro-differential  equations,  with  constant  coef- 
ficients. When  inputs  to  this  system  are  mentioned,  the 
terms  "step",  "ramp",  and  "parabola"  of  input  will  be  used, 
A  step  of  the  input  quantity  is  that  input  which  yields  zero 
error  when  there  is  unity  feedback  from  output  to  input.   For 
a  Type  I  positional  servo,  a  step  would  be  a  step  of  position. 
For  a  Type  II  velocity  servo,  the  step  would  be  a  step  of 
velocity,  etc.   By  treating  such  a  general  system,  the  ensuing 
results  can  be  applied  to  nearly  all  practical  systems. 

This  analytical  compensation  design  technique  uses  as  its 
cornerstone  the  characteristic  equation  of  the  system.   After 
a  brief  discussion  of  the  methods  for  choosing  direct  path 


components,  a  method  for  forming  this  characteristic 
equation  is  described.   The  system  is  represented  in  a 
standard  block  diagramo   Using  circuit  analysis  techniques, 
the  interaction  of  the  blocks,  path  and  loops  in  this  dia- 
gram is  shown  in  a  characteristic  determinant  of  the  system. 
It  is  from  this  determinant  that  the  characteristic  equation 
is  written. 

The  next  step  in  the  procedure  relates  performance  speci- 
fications to  the  coefficients  of  the  characteristic  equation. 
Static  specifications  such  as  "following  error"  and  output 
in  response  to  a  load  torque,  place  limits  on  coefficients 
or  on  ratios  of  coefficients  in  the  characteristic  equation. 
Because  of  the  vast  amount  of  work  which  has  been  done  on 
simple  second  order  systems,  most  dynamic  specifications 
are  expressible  in  terms  of  the  parameters  of  such  systems. 
The  compensated  system  is  made  to  have  a  response  not  sig- 
nificantly different  from  the  specified  second  order  response. 
This  is  done  by  making  the  quadratic  characteristic  equation 
of  the  second  order  system  a  factor  of  a  model  characteristic 
equation  for  the  compensated  system.   All  other  roots  of  this 
model  equation  are  forced  to  be  so  large  that  the  second  order 
roots  dominate  the  response. 

The  effect  of  all  possible  compensators  on  the  coeffi- 
cients of  the  characteristic  equation  is  found  through  the 
use  of  the  characteristic  determinant.  All  possible  paths  in 
the  system  are  uniquely  represented  in  this  determinant. 
Closing  any  of  these  paths  with  a  compensator  adds  a  cofactor 
to  the  existing  determinant,   A  rapid  method  for  determining 


the  effect  of  this  path  closing  on  the  coefficients  of  the 
characteristic  equation  is  presented.  This  information  is 
displayed  in  a  Location  Matrix. 

At  this  stage  of  the  design  process,  the  designer  has  a 
model  characteristic  equation  representing  the  specifications, 
and  an  actual  characteristic  equation  representing  the  uncom- 
pensated system.   A  comparison  of  these  shows  which  coefficients 
must  be  altered  by  compensation.   The  Location  Matrix  indicates 
what  compensator,  in  what  location,  will  effect  the  required 
change . 

This  thesis  confines  itself  to  the  problem  of  assuring 
dominance  of  a  desired  second  order  response.   The  technique 
itself  is  not  so  limited,  and  could  be  used  to  force  a  system 
to  perform  like  any  desired  third  or  higher  order  systemo 


CHAPTER  2 
BASIC  DESIGN  CONSIDERATIONS 

lo   Introduction 

Some  servoraechanisms  perform  a  computing  function. 
Examples  of  this  type  of  servo  are  found  in  fire  control 
systems.   In  these  servos,  the  output  or  outputs  may  be  com- 
plicated mathematical  functions  of  several  inputs.   The  mathe- 
matical relationships  to  be  instrumented  will  dictate  the 
choice  of  equipment  to  be  used.   The  type  of  servo  of  primary 
interest  in  this  thesis,  however,  is  one  which  drives  a  load 
so  that  its  action  duplicates  the  action  represented  by  the 
input  signal.   An  example  of  thia.  type  of  servo  is  th®  drive 
system  for  a  gun.   Its  input  is  a  voltage  or  shaft  rotation 
which  represents  the  desired  orientation  of  the  gun  barrel. 
It  is  the  function  of  the  servo  to  control  the  power  that 
drives  the  gun  mount  so  that  the  desired  gun  barrel  orien- 
tation is  achieved.   The  preliminary  choice  of  equipment  to 
perform  this  type  of  function  is  considered  in  this  chapter. 
These  considerations  result  in  a  basic  system  with  sufficient 
power  and  accuracy  to  meet  the  fundamental  requirements  3e% 
by  the  specifications.   Subsequent  chapters  deal  with  modi- 
fying the  basic  design  to  meet  specifications  concerning  the 
more  detailed  nature  of  the  servo's  response. 


2.   Components  of  a  Typical  Servomechanlsm 

Figure  II-l  is  a  diagram  of  a  typical  servomechanism.   Its 
purpose  might  be  to  control  the  angular  position  and/or  velocity 
of  a  shaft.   The  shaft  is  represented  by  the  block  labeled 
"Load".   The  input  or  command  might  be  the  angular  position 
of  a  shaft  from  a  computer.   The  command  receiver  and  the  out- 
put measuring  device  convert  the  input  and  output  into  quanti- 
ties which  can  be  easily  compared.   The  comparator  subtracts 
these  quantities  and  produces  a  signal  proportional  to  the  dis- 
parity between  input  and  output;  i.e.,  error.   Since  this  error 
may  be  a  small  quantity,  it  must  be  amplified.   It  is  then' 
applied  to  a  controller  which  actuates  the  prime  mover.   The 
prime  mover  drives  the  load  through  the  gear  train  so  as  to 
null  the  error. 


Input  >  Command_jpomp 
iRec'  v'ri  oatQi 


Error 


jAmp  li4_jCont-_JPrime 
fier        roller    piovei 


Output 

Measuring 

Device 


Load 


Output 


Figure    II-l 

A   servomechanism  can  be   considered   to    consist   of  paths 
through  which  information  flgws.      Those  paths   which  conduct 
information  in  the    direction  from  input    to   output   are   termed 
forward  paths.      Those  which   conduct   in   the   opposite   direction 
are    termed  reverse   paths.      In  essence,    the  forward  paths    of   a 
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servo  transform  the  in^ut  into  the  output.   The  reverse,  or 
feedback,  paths  insure  the  proper  degree  of  correspondence 
between  input  and  output. 

In  a  relatively  simple  system,  there  is  usually  one  for- 
ward path.  However,  other  forward  paths  may  result  from  com- 
pensation considerations.  To  avoid  ambiguity,  the  forward 
path  which  performs  the  primary  function  of  changing  input  t^ 
output  will  be  designated  the  main  path.  Other  forward  paths 
will  be  designated  feed-forward  paths. 

3 •   Some  Considerations  Governing  Choice  of  Components 

The  need  for  a  servo  implies  that  the  element  which  generates 
the  input  has  insufficient  power  to  control  the  load.   Since 
the  purpose  of  the  servo  is  to  supply  the  needed  power,  choice 
of  the  prime  mover  is  a  logical  starting  point  in  selecting 
main  path  components.   The  characteristics  of  the  load  will 
be  known.   These  characteristics  may  be  expressed  in  terms  of 
its  inertia,  the  viscous  friction  to  which  it  will  be  subjected, 
external  torques  which  may  act  on  it  and  the  like.   The  speci- 
fications will  also  state  the  velocities  and  accelerations  re- 
quired of  the  load.   This  information  will  set  a  lower  limit 
on  the  power  and  torque  capacity  required  of  the  prime  mover. 
The  power  supply  available  will  further  limit  the  choice.   Other 
requirements  on  the  prime  mover  will  arise  from  cost,  weight 
and  size  limitations,  life  expectancy,  etc.   With  these  limits 
in  mind,  the  designer  must  then  select  a  prime  mover  from  those 
commercially  available  or  else  design  his  own. 

Choice  of  the  prime  mover  will  normally  place  limits  within 


which  a  gear  ratio  must  be  chosen.   Ahrendt   discusses  the 
various  factors  which  influence  this  choice.   Such  factors 
as  low  cost  or  smoothness  of  output  at  low  velocities  may 
place  conflicting  demands  on  the  choice  of  gear  ratio.   Also, 
the  inpedance  matching  aspect  of  gear  ratio  selection  may  be 
of  great  importance.   In  order  to  choose  the  optimum  gear  ratio, 
the  designer  must  have  a  thorough-going  knowledge  of  the  speci- 
fications as  well  as  of  the  characteristics  of  the  prime  mover. 

Choice  of  the  prime  mover  will  probably  narrowly  limit 
choice  of  the  controller.   If  a  variety  of  controllers  is 
available,  the  considerations  of  fineness  of  control  versus 
cost,  weight,  etc,  must  be  brought  to  bear. 

Insufficient  information  will  be  available  in  the  speci- 
fications to  allow  choice  of  the  amplifier  to  be  used.   The 
gain  to  be  provided  will  be  determined  by  compensation  considera- 
tions.  This  will  not  be  known  until  later  in  the  design  process. 
When  it  has  been  determined,  the  amplifier  may  be  chosen.   In 
addition  to  the  conditions  imposed  on  amplifier  selection  by 
such  factors  as  available  power  supply,  noise  generation,  and 
operating  environment,  the  designer  must  bear  in  mind  the  impor- 
tance of  linearity.   Although  sufficient  band  width  is  not  diffi- 
cult to  obtain,  saturation  effects  may  present  a  problem.   It 
is  uneconomical  for  a  system' s  amplifier  to  saturate  before  its 
prime  mover  does,  as  this  means  unusable  power  has  been  purchased, 

The  choice  of  comparator,  measuring  device  and  command 
receiver  will  be  a  question  of  obtaining  the  best  linearity  and 
accuracy  within  the  limits  set  by  cost,  weight  and  size  speci- 
fications.  These  devices  are  described  in  detail  in  several 
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12  3 
servome Chan ism  texts.  *  ' 

i^.   Determination  of  Component  Transfer  Functions 

Having  chosen  the  main  path  components,  their  transfer 
functions  must  then  be  determined.   For  components  which  have 
their  output  independent  of  the  frequency  of  the  input,  this 
merely  involves  determining  their  sensitivity.   It  is  desirable 
that  this  sensitivity  will  be  constant  over  the  range  of  interest 
in  the  problem.   If  sensitivity  is  not  constant,  it  must  be 
assumed  so  in  the  remainder  of  the  problem  by  assigning  it  a 
representative  value. 

The  transfer  functions  of  frequency  dependent  components 
may  be  determined  mathematically  by  application  of  basic  a-c 
circuit  theory  and  the  Laplace  Transformation.   This  technique 
is  described  in  detail  by  Thaler^,   If  the  electrical  parameters 
of  the  component  are  unknown,  it  may  be  necessary  to  conduct 
tests  on  the  equipment  to  determine  its  transfer  function. 
These  may  take  the  form  of  frequency  or  transient  response  tests. 
The  results  of  these  tests  will  permit  derivation  of  the  trans- 
fer function.   Examples  of  this  method  are  presented  in  refer- 
ences 5  and  6. 

Having  obtained  the  transfer  functions  of  the  equipment 
to  be  used,  it  is  possible  to  express  the  system  in  this  stage 
of  its  design  by  a  signal  flow  chart  or  standard  block  diagram. 
This  phase  of  design  is  discussed  in  the  next  chapter. 
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CHAPTFR  3 
FORMATION  OF  THE  CHARACTERISTIC  EQUATION 

3.1  Introduction.   The  characteristic  equation  is  an  equation 
in  s  domain  which  represents  the  integro-differential  equa- 
tions in  time  space  which  describe  the  flow  of  information 

7 
through  a  system.    The  roots  of  this  equation  allow  the 

determination  of  the  transient  performance  and  the  stability 

of  the  physical  system  so  described.   The  purpose  of  this 

chapter  is  to  develop  a  standard  method  of  forming  this 

characteristic  equation  for  a  multiple-loop  servo  system. 

3.2  Standard  Block  Diagram.   A  set  of  transfer  functions  are 
developed,  as  described  in  Chapter  2,  for  elements  of  the 
servomechanism.   The  block  diagram  portrays  the  interconnec- 
tion of  these  various  elements  to  form  the  overall  system. 

The  very  nature  of  the  block  diagram  assumes  that  the  transfer 
function  of  the  tandem  conbination  of  two  blocks  is  the  product 

o 

of  the  individual  transfer  functions.    The  overall  perfor- 
mance of  the  system  is  found  by  combining  these  blocks 

9  10 

according  to  the  standard  rules  for  their  manipulation.  ' 

Unfortunately  for  the  designers,  the  combining  of  the  blocks 
destroys  their  individuality  and  the  identity  of  the  mechanism 
or  circuit  they  represent »  What  is  needed  is  an  alternate 
representation  of  the  system  which  pictures  it  in  more  detail 
than  the  reduced  (combined)  block  diagram.   This  could  take 
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a  new  form,  as  in  the  signal-flow  diagrams  originated  by 
S.  J.  Mason,    or  it  could  be  a  more  ordered  manipulation  of 
the  block  diagram,  as  proposed  by  Yaohan  Chu.    It  is  to 
the  latter  method  that  we  now  devote  our  attention. 

Certain  advantages  accrue  from  standardizing  the  block 
diagram.   In  this  regard  the  following  symbols  are  defined: 

a.  The  Nodal  Point  or  Node:   These  are  summation  points 
in  the  main  path  of  the  system.   Each  node  may  receive  a 
number  of  signals,  but  delivers  only  one.   The  nodes  will  be 
lettered  in  sequence,  a,  b,  c,  .  .  n,  from  left  (input)  to 
right  (output). 

b.  The  Pick-off  Point:   This  is  a  point  at  which  there 
is  one  incoming  signal  and  two  or  more  outgoing  signals.   The 
outgoing  signals  are  identical  with  the  incoming  one. 

c.  The  Block  Transfer  Function;   The  transfer  functions 
in  the  blocks  will  be  designated  by  double  subscript  notation. 
This  notation  is  defined  for  the  three  types  of  blocks  as 
follows: 

1)  Direct  Path  Block:   These  blocks  are  noted  G  „  where 
the  "n"  refers  to  the  node  preceding  the  block.   "N"  is  the 
upper  case  form  of  the  same  letter. 

2)  Feedback  Block:   The  feedback  blocks  are  subscripted 
GjT  .   The  "N"  is  the  same  as  the  upper  case  subscript  of  the 
direct  path  block  from  which  they  receive  their  input.   The 
"m"  is  the  designation  of  the  node  to  which  they  feed. 

3)  Feed-forward  Block:   Such  blocks  contain  the  symbol 
G  „.   The  lower  case  subscript  refers  to  the  main  path  node 
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preceding  the  block.   The  upper  case  subscript  is  the  same  as 

the  upper  case  subscript  on  the  main  path  block  which  feeds 

into  the  same  node . 

As  examples  of  the  above,  consider  the  following  three  node 

system: 

Figure  3-1 


Output 


@,  @>  and  (c)  are  nodes. 


h. 


symbolizes  a  pickoff  point 


Although  the  system  above  is  represented  as  having  but  one 
input  and  one  output,  there  is  no  limit  to  the  number  of 
either.   Internal  noise  generated  in  one  block  may  be  considered 
as  an  external  input  to  the  node  following  that  block.   If 
there  is  more  than  one  input,  each  is  designated  by  a  Roman 
numeral.   If  there  is  a  transfer  function  between  an  input  and 
its  entering  node  it  will  be  subscripted  as  follows; 


(Input  "I") ) 


la 


Ih 


3.3  Mathematical  Analysis  Based  on  the  Standard  Block  Diagram. 

3.3.1  Formation  of  the  Characteristic  Determinant:   We 
can  obtain  a  set  of  simultaneous  equations  for  the  system 
represented  by  a  standard  block  diagram  by  applying  the  network 

1^5 

analysis  technique  of  Guillemin.  -^   This  entails  writing  the 
equation  for  the  signal  summation  at  each  node.   If  we  iden- 
tify the  signal  out  of  a  node  with  the  lower  case  letter  of 
the  node,  we  have  the  following  equations  in  our  examples 
At  Node  a: 

«  =  I(°la)  -  =<°ccGca>  -  ^(GaA°Aa'  '  ^.(°bB°Ba) 

At  Node  b: 
b  =  a(G^^)  -  bCGj^gGgj,)  -  o(G^(,°Cb^ 

And,  at  Node  c: 
c  =  a(G^3)  +  bCGj^g)  -  o(G^(,G(,^) 

Rearranging  these  equations  we  haves 

(1  *   °aA°Aa'^  ^  ^SbSJ"  *    (°cC°Ca^=  .=  ^^la^^ 

(-G^^)a  +  (1  +  G^3G3,^)b  h-  (G^cO^^^)c  =   0 

(-G^B)a  +  (-G^£)b  +  (1  +  G^c°Cc'=  =  ° 

Inspection  of  these  equations  reveals  a  symmetry  of  arrarlge- 
ment  which  would  permit  writing  the  ecjuations  almost  without 
reference  to  the  block  diagram.   All  the  forcing  functions  are 
on  the  right.   On  the  left  are  terms  which  contain  components 
and  signals  within  the  system.  We  can  express  these  equations 
in  matrix  notation  as  follows; 
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^  ■"  ^aA^Aa      ^bB^Ba     ^cC^Ca 


-G 


aA 


^  ■"  ^bB^Bb  ^cC^Cc 


-G 


aB 


-G 


bB 


1  ^  ^C^Cc 


^(^la^ 


We  call  the  determinant  of  the  matrix  of  the  coefficients 
the  characteristic  determinant  of  the  system.   To  this 
determinant  we  assign  the  symbol  A. 

3. 3*2  Manipulations  of  the  Characteristic  Determinant* 
Independent  of  its  physical  significance,  the  determinant  of 
the  coefficients  of  a  set  of  linear  equations  has  certain 
properties.   These  properties  are  illustrated  in  the  tech- 
nique of  equation  solving  with  determinants  and  in  the  rules 
for  expansion  of  a  determinant. 

The  solution  of  the  equations  is  accomplished  as  follows: 

1)  Replace  the  column  corresponding  to  the  desired 
unknown  with  the  column  of  forcing  functions.   The  resulting 
determinant,  divided  by  /\,    is  the  value  of  that  unknown  in 
response  to  these  inputs. 

2)  Since  the  determinant  represents  a  linear  system,  the 
principle  of  superposition  is  valid.  Hence,  in  the  case  of 
multiple  inputs,  the  output  resulting  from  each  may  be  com- 
puted separately  and  the  individual  results  added. 

3)  The  output  obtained  by  replacing  the  last  column  of 
the  determinant  with  the  input  column  is  the  signal  from  the 
last  node.   It  must  be  multiplied  by  the  performance  function 
following  the  last  node  to  obtain  the  system  output. 


16 

The  characteristic  (homogeneous)  equation  describing  the 
system  is  obtained  by  setting  the  characteristic  determinant 
equal  to  zero. 

The  normal  rules  for  expanding  a  determinant  apply  to  the 
characteristic  determinant.   In  general  they  require  the  expan- 
sion by  co-factors  of  some  row  or  column.   Since  we  will 
usually  be  interested  in  the  final  output,  we  should  expand 
by  elements  of  the  last  column.   To  the  co-factor  of  the  ith 
row  and  jth  column  we  assign  the  symbol  /!-«■!•   Similarly  the 
element  in  the  ith  row  and  jth  column  is  called  G .  . .  Using 
this  notation,  the  expansion  of  an  "n  by  n"  characteristic 
determinant  by  elements  of  the  last  (nth)  column  is: 

n  .    .       . 

As  an  example  of  the  procedure  described  above  for  finding 
system  output  in  response  to  a  given  input  consider  the 
following: 

Input  to  node  "a"  is   I,  G^^  =  1 


«cc)I 


Output  = 


■°aA    1  *  SfiOBb 


•°aB :Sb 


H 


°"*P"*  =  ^(°eC^^°aA°bB  ^  °aB  ""   °aA°bB°Bb7 

3.30  Properties  of  the  Characteristic  Determinant.   In 
addition  to  properties  inherent  in  all  determinants  the  charac' 
teristic  determinant  has  featui-es  peculiar  to  the  system  it 
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represents. 

1)  The  characteristic  equation,  obtained  by  equating  /\ 
to  zero,  is  the  sum  of  unity  plus  a  number  of  closed  loop 
performance  functions.   By  the  nature  of  the  elements  which 
comprise  the  determinant,  these  functions  will  be  ratios  of 
polynomials  in  s.   When  this  sum  is  placed  over  a  common 
denominator,  its  numerator  will  be  the  usable  form  of  the 
characteristic  equation  of  the  system. 

2)  Each  element  of  the  characteristic  determinant 
corresponds  to  a  path  in  the  standard  block  diagram. 

3)  The  diagonal  elements  of  the  characteristic  deter- 
minant correspond  to  the  paths  of  inner  loops  of  the  nodes. 
The  elements  above  these  diagonal  elements  are  all  feedback 
paths.  Those  below  the  diagonal  are  feed-forward  or  direct 
paths. 

h)  Each  feedback  and  feed-forward  performance  function 
appears  in  only  one  element  of  the  determinant.  Each  direct 
path  performance  function  may  appear  in  more  than  one  element. 

We  now  have  at  our  disposal  a  standard  block  diagram. 
This  block  diagram  maintains  the  integrity  of  each  element 
comprising  the  system  it  represents.  We  form  this  block 
diagram  to  aid  us  in  compensating  the  system.   To  this  end 
we  need  only  include  those  nodes  or  pick-offs  which  represent 
points  in  the  system  physically  usable  for  compensation.   All 
performance  functions  between  such  points  should  be  combined 
into  one  block  to  simplify  the  ensuing  algebra.   The  expansion 
of  the  standard  determinant  representing  the  block  diagram  gives 
the  characteristic  equation  of  the  system.   Chapter  h   will  relate 
performance  specifications  to  this  characteristic  equation. 
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CHAPTER  h 

THE  EFFECT  OF  PERFORMANCE  SPECIFICATIONS  ON  THE 
COEFFICIENTS  OF  TIffi  CHARACTERISTIC  EQUATION 

1.  Introduction 

The  performance  specifications  state  the  desired  system 
output  in  response  to  actuating  or  disturbing  inputs.   The 
output  is  specified  in  terms  of  its  steady  state  value  and  its 
transient  behavior. 

The  characteristic  equation  which  was  evolved  in  the  pre- 
vious chapter  describes  the  system  performance  in  the  s  domain. 
The  standard  analysis  approach  is  to  find  the  roots  of  this 
characteristic  equation  and  transform  it  to  the  time  domain 
where  actual  and  specified  performance  can  be  compared. 

The  purpose  of  this  chapter  is  to  relate  the*  specifications 
directly  to  the  coefficients  of  the  characteristic  equation. 
This  eliminates  the  need  for  finding  the  roots,  plotting  the 
response,  and  other  laborious  analysis  techniques. 

2.  Static  Performance  Specifications 

The  static  performance  specifications  describe  the  steady 
state  value  of  the  system  output »   In  response  to  an  actuating 
input  the  specification  will  normally  state  an  allowable  error. 
In  response  to  a  step  we  would  expect  no  error.  A  following 
error  will  be  stated  for  the  response  to  a  ramp  input.   A 
maximum  allowable  output  is  usually  stated  for  the  response  to 
a  disturbing  input.   These  disturbances  might  be  noise  or  some 
form  of  load  torque. 
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Evaluation  of  the  steady  state  output  to  a  step,  ramp  or 
parabola  of  input  position  is  most  useful  in  studying  the  effect 
of  disturbance  signals.   As  established  in  Chu' s  dissertation 
the  performance  function  relating  the  input  at  any  node  with 
the  output  at  any  other  node  is 

Output   ^  ^JjAiJ 

Input        A  ^^   ' 

where:      i^  node  into  which  input  is  fed, 

j^  node  just  preceding  output  block, 
G.J  ^    performance  function  of  output  block, 
/\ .  .  ■^  cofactor  of  the  i   row  and  the  j    column 
of  the  standard  determinant. 
The  output  in  response  to  a  given  input,  I(s)  is 
Ks)  G   An 

A 

The  steady  state  value  of  this  output  is  obtained  by  applying 
the  final  value  theorem  as  follows: 


Output 

"    s  -  0 


i'^)  °.ijAi.i' 


A 


(1^.3 


o.      A  _   Characteristic  Equation      /i  is 
Since  l\  -     g^— ^ {kA) 

where  CD  a^  common  denominator  of  the  ratios  Qf 
polynomials  making  up  the   expansion  of  /\. 


■    Characteristic  Eqn.J 


0utput^3  =  Limit  I  nv.o^..^i::.-otj.  ^L^    I      (^-5) 


^0 


For   a    step    of   disturbance,    i.e.      I(s)   =  a/s 


Output        =  Limit 
^^        3   ->  0 


o 


since  a 
o 


is  the  limit  of  the  characteristic  equation  as  s  goes  to  zero. 
For  a  specified  maximum  output  over  input  ratio,  k,  where 

IT 

k  >  Limit   [gjjAij  CD  a]  ik-.e) 

A  close  examination  of  the  equation  above  will  reveal  the 

effects  of  such  a  specification  on  the  a   coefficient  of  the 

o 

characteristic  equation. 

Case  I  -  If  k  is  a  number  not  equal  to  zero,  then  a  is 
specified  as  being  greater  than  the  Limi 


t    KuAlj  cd/\1 


This  may  or  may  not  limit  the  size  of  a   since  this  limit  may 

be  identically  zero. 

Case  II  -  If  k  is  specified  as  zero  then  the 

must  be  zero  since  a^  cannot  be 
o 


lite . 

In  Case  I  where  the  Limit  |  G  .  t/\  .  •  CD  a    i 
s  -^  0  L  J^   ^J   ^-AJ 


infinite 

is  a  constant 


we  have  set  one  condition  on  a  which  may  be  expressed  in 
equation  form  as: 

a   >  N 
Note  that  the  Limit     G^t  Z\-  •  ^^A  ^^7   ^®  constant  containing 

variable  gains  of  the  uncompensated, system.   The  condition  on 
a  may  be  met  by  setting  these  gains  which  determine  the  stiff- 
ness of  the  system.   This  will  then  place  a  lower  limit  on 
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these  gains  and  restrict  their  use  in  future  compensation. 

In  Case  II,  if  the  Limit   G .  t  A  4  4  CDa   ^   0>   then  an 
s  ->  0   L  ^"^   ^J    ^J 

s  must  be  introduced  into  the  numerator  of  /\. .  which  does  not 

cancel  an  s  in  its  denominator;  or  an  s  must  be  placed  in  CD/\ 

which  is  not  also  placed  in  the  denominator  of  /\,  .,        This 
will  raise  the  order  of  the  characteristic  equation.   The  former 
may  be  accomplished  by  introducing  derivative  feedback  from 
any  block  before  the  output  block,  to  the  actuating  input's 
node.   The  introduction  of  the  excess  s  in  CDa   is  brought 
about  by  putting  integration  feedback  from  other  than  the  out- 
put block  to  the  input  node  to  which  the  disturbing  inputs 
are  applied. 

To  illustrate  the  reasoning  behind  the  statements  above 
consider  a  four  node  system  represented  by  the  following  stan- 
dard block  diagram: 


'~m 


aA 


-<§>- 


'bB 


^ccj- 


->(S>— 


•dD 


DA 


The  standard  determinant  for  this  system  is: 


A 


^aA 
0    -G 


bB 
0    -G 


cC 


^dD^Da 
0 

0 

1 


For  a  disturbing  input  at  node 


we  wish  the  output  to 


be  zero;  but  have  postulated  that   Limit   G,_^  ZX"^)!  ^A    ^   ^' 
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We  need  to  introduce  an  s  into  one  of  the  terms  of  this  express- 
ion which  will  drive  the  response  to  zero  without  affecting  the 
response  to  the  actuating  input. 

The  steady  state  value  of  the  response  to  an  actuating 
step  input,  A/s,  is  (by  equation  1^,5)' 


1^   [GdiAi^ci^A] 


Limit   A   iG.^Ani.  CD  a  I  =  A 
s  -♦  0 


Since  G,_.  and  CD  a   are  common  to  both  expressions  a  change  in 
either  will  affect  both  responses.   If  the  s  is  introduced  in 
CD/\   it  must  also  appear  in  the  denominator  of  /\-.  i   so  that 
its  effect  is  not  felt  in  the  case  of  the  actuating  input.   It 
must  not,  however,  appear  in  the  denominator  of  </\^i.«   The  2nd 
and  I|.th  rows  of  /\   are  contained  in  both  cofactors  so  they 
cannot  be  used.   The  3rd  row  appears  in  /\-.  i  but  not  in  /\-:^\   . 
The  only  locations  usable  for  an  integration  feedback  which 
would  have  the  desired  result  are  positions  in  the  3rd  row. 
In  our  example  the  only  such  feedback  position  is  around  G  „. 

The  other  method  of  introducing  an  s  in  the  expression  for 
steady  state  output  in  response  to  a  disturbance,  without  alter- 
ing the  desired  actuated  response,  is  to  put  an  s  in  the  numer- 
ator of  /\ri\    which  does  not  appear  in  the  numerator  of  /\->\   • 
This  is  accomplished  by  inserting  derivative  feedback  in  a 
first  row  position.   The  1st  row  is  the  only  one  in  /\-:,\    which 
is  not  also  in  /\-,\   . 

In  higher  order  systems  there  may  be  peculiarities  of  the 
determinant  which  would  allow  other  positions  to  be  used,  or 
which  would  eliminate  some  of  the  above  positions.   Those 
explained  above  are  the  only  positions  possible  if  all  terms 
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Ip.   the  determinant  are  present. 

Computation  of  the  output  quantity  is  of  prime  import- 
ance in  devising  methods  of  making  a  system  insensitive  to 
disturbing  inputs.   In  response  to  an  actuating  input  we 
might  better  express  steady  state  performance  in  terms  of 
error  since  this  permits  an  immediate  and  direct  check  on 
the  accuracy  of  the  system.   If  error  is  defined  as  the 
difference  between  the  actuating  input  and  its  desired  out- 
put, then  the  error  equation  is: 

E  =  Input  -  Output  (^.7) 

This  may  be  rewritten  as: 

Input.         (^.8) 


Output 
Input 


Expressing  Output/input  as  per  equationC+.l)  we  have: 


E 


mAj 


l(s)  (^.9) 


A 

By  the  final  value  theorem,  and  again  expressing  /\     as  in 
equation  (^.^),  we  have 

Error^g  =  Limit   sl(s)  Ichar.Eqn.  -  G .  jAijCDA  I   (^.10) 


Char . 


^  -*  ^  '        Char.  Eqn. 


In  response  to  I(s)  =  A/ss 


.mit  A   Char  Eqn  -  O.jA.. 

^   0      hi — 12. 

L      Char,  Eqn. 


Error^g  =  Limit  A   Char  Eqn  -  G-^jAj-j  ^'^/\  I  (^-ID 


Since  the  steady  state  error  in  response  to  a  step  will  be 
defined  as  zero  we  can  write: 


2h 


Limit  A  Char.Fqn.  -  G-jA^.  CD  a 

0  =   ^  "^  ^ == -^  (^-.12) 

^  Limit    Char.  Fqn.  ^^.-l^; 

s  -^-  0 
Since  Limit  Characteristic  Equation!   =  a^  0,  then: 

^°  ^  s'fo  hj^iJ  ™a1       ^^-"^ 

This  defines  another  condition  on  a  ,  but  one  which  will 
automatically  be  met  if  there  is  unity  feedback  from  the 
output  block  to  the  actuating  input  node. 

In  response  to  an  input  ramp,  i.e.  I(s)  =  B/s  : 


Error   =  Limit  B 
^^   s  ->  0  s 


Char.  Eqn.  -  G.j/\j_-j  CD  a 
Char.  Eqn. 

(^.1^) 


For  zero  error  in  response  to  a  step,  the  Limit  of  the 

s  -*  0 

Characteristic  Equation  equals  the  Limit  G.j/\_.  .  CD  a  =  a 

Hence  the  constant  parts  of  the  numerator  of  equation  C+-.l^) 
are  equal.   This  means  that  s  is  a  factor  of  the  numerator. 


t  1        fais   ^  a^      -     G.jAj,   CD^^j   {h.. 


Error^g  =   Limit   B        fai  s   +  a^      -     G.^A^.   CD  a  I    (^.15) 


Error  _   ai    .. 

B     ^^     "   a"  s  -;   u 

o 

This  equation  puts  a  limit  on  the  ratio  of  a-j/a^  in  the 
compensated  characteristic  equation. 

This  ratio  m^|r  not  be  compatible  with  the  a-ji/a^  ratios 
of  any  of  the  model  equations  formed  from  dynamic  specifi- 
cations.  The  restriction  on  a-j_/a^  may  be  relieved  by  makiftg 
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the  ramp  following  error  identically  zero,  i.e.,  raise 
system  to  Type  II.    ' 

All  the  conditions  on  equation  order,  compensator  loca- 
tion and  coefficient  magnitude  which  arise  from  consideration 
of  static  specifications  should  be  held  for  simultaneous 
solution  with  similar  conditions  yielded  by  the  dynamics  of 
the  problem. 

C.   Dynamic  Performance  Specifications 

The  dynamic  performance  specifications  describe  the 
desired  transient  performance  of  the  system.  Since  they  are 
meaningless  in  an  unstable  system,  they  all  imply  stability 
even  if  they  do  not  state  it. 

As  opposed  to  the  static  specifications,  where  there 
is  some  uniformity,  there  are  nearly  as  many  ways  of  express- 
ing a  desired  transient  behavior  as  there  are  Process  Engin- 

12 
eers.   The  following  summary   ,  while  not  intended  to  be  all 

inclusive,  lists  the  most  common  specifications  and  discusses 

how  to  calculate  them. 

THE  TEN  MOST  COMMON  DYNAMIC  PERFORMANCE  SPECIFICATIONS 


NAME 


TYPE 


MFTHOD  OF  COMPUTATION 


1 

Gain 
Margin 


Stability  specifi 
cation 


Compute  maximum  gain 
for  stability  by  Routh- 
Hurwitz  criterion.  Ratio 
of  maximum  stable  gain 
to  actual  gain  is  gain 
margin. 
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2 
Phase 

Margin 


Frequency  domain 
stability  specifi- 
cation. 


3 

M  Peak 


Frequency  domain 
stability  spec- 
ification. 


On  Nyquist  diagram,  the 
angle  between  curve  and 
negative  real  axis  at  the 
unit  circle.   On  Bode 
plot,  find  00^  (frequency 
of  unity  gain),  compute 
phase  shift  at  co^,  and 
subtract  from  l80  deg. 
For  a  well  behaved  system, 
a  phase  margin  of  ^5  deg. 
means  one  overshoot  in 
response  to  a  step. 

The  M  curve  is  the  mag- 
nitude of  the  closed- 
loop  response  vs. 
frequency.   It  can  be 
calculated  from  open-loop 
transfer  function  or  from 
values  picked  off  Nyquist 
plot.  M  peak  is  defined 
as  the  peak  of  closed- 
loop  frequency  response 
normalized  to  low  fre- 
quency value. 

M  =     A 

1  -  A 

Circles  of  constant  M 

can  be  drawn  on  Nyquist 


Damping 
Ratio,  C 


Stability  speC' 
if ication. 
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plot.   The  locus  of  the 
system  must  remain  out- 
side given  M  peak  circle. 
(An  M  peak  of  1.^  usually 
means  one  overshoot) 

Can  be  calculated  from 
closed-loop  transfer 
function  or  found  directly 
from  root-locus  plot. 
Damping  ratio  is  defined 
as  C  in  a  quadratic  term 

of  the  form: 

2  2 

s  +  ^^^rf  "^  ^n 

On  the  root-locus  dia- 
gram, C  equals  cos  9 
where  9  is  the  angle 
between  the  horizontal 
and  the  vector  directed 
toward  the  gain  point 
on  the  locus. 


5 

Damping 
Factor,  Co) 


Stability 
Specification 


n 


or  Decrement 
Factor 


May  be  found  directly 
from  the  root-locus. 
Damping  factor  is  defined 
by  the  factors  of  a  quad- 
ratic system,  as  was 
damping  ratio 5  it  is  the 
real  portion  of  the  roots. 
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6  Direct  evaluation 

Percent     of  system 
Overshoot   stability 


7  Frequency  domain 

Bandwidth   speed-of -response 

specification 


Where  the  imaginary  por- 
tion of  the  root  deter- 
mines the  frequency  of 
the  damped  sinusoidal 
response  of  a  system, 
the  real  portion  deter- 
mines the  rate  of  decay. 

May  be  calculated  from 
response  of  system  to 
step  function  input. 
Defined  as  peak  of  the 
response  to  a  step-func- 
tion input  as  a  percent- 
age of  final  value o 

Usually  defined  as  the 
frequency  at  which  the 
closed-loop  frequency 
response  falls  to  l//2 
its  low  (zero)  frequency 
value.   Bandwidth  is 
available  directly  from 
the  M  curve  of  the  system 
An  analogous  definition 
would  be  the  crossover 
frequency  of  the  system. 
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8  Speed-of -response 

Rise  Time   specification 


One  of  the  simplest 
definitions  of  rise  time 
is  l/bandwidth,  another 

^^  ■'■/"^  peak*  ^^®®  ^^°^® 
is  also  defined  as  time 

to  first  zero  of  error. 


9  Overall  performance 

"Settle-Out"  Specification 

Time  or 

"Synchronizing"' 

Time 


Settle -out  time  can  be 
calculated  analytically 
by  the  inverse  Laplace 
transform,  or  by  picking 
values  off  the  root  locus 
plot.  This  specification 
is  most  convenient  if  an 
analog  computer  is  avail- 
able. Settle-out  time 
is  defined  as  the  time 
required  for  system 
response  to  a  step-func- 
tion input  to  approach 
and  remain  within  a 
given  tolerance  of  the 
final  value.   It  is 
usually  necessary  to 
specify  maximum  over- 
shoot and  steady-state 
error  in  addition  to 
settle-out  time. 
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10 

ITAE  or  Integrated 

Value  of  the 

Product  of  Time 

and  Absolute 

Value 


Overall 

Performance 

Specification 


The  ITAE  specification 

is  defined  by  the  integral 

CO 

J   tietdt 
0 

Minimize  the  value  of 

this  integral  for  optimum 

system  performance.  This 

specification  yields  a 

number  that  depends  on 

system  parameters.   It 

would  probably  be  an 

improvement  to  generalize 

it  by  normalizing  to  the 

rated  output  quantity. 

The  dimensions  would  then 

2 
be  sec  .  Another  poss- 
ible generalization 
would  be  to  integrate 
to  some  arbitrary  large 
number,  such  as  10  time 
constants,  rather  than 
to  infinity.  This  would 
permit  the  specification 
to  be  applied  to  responses  for 
which  the  system  yields 
a  finite  final  error. 
ITAE  is  one  of  several 
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attempts  to  define  an 
overall  figure  of  merit 
for  system  operation. 
Rather  than  simply  sum- 
ming the  total  error  of 
a  system,  the  error  is 
progressively  weighted 
more  heavily  as  time  goes 
on.  This  puts  a  premium 
on  rapid,  accurate  settle- 
out,  and  allows  for 
unavoidable  initial 
large  error. 


Due  to  the  complexity  and  variety  of  responses  of  sys- 
tems of  order  higher  than  second,  very  little  information 
is  available  in  the  literature  concerning  the  characteristics 
of  their  response.   However,  the  response  of  a  second  order 
system  is  easily  computed,  fairly  simple,  and  can  be  charac- 
terized by  two  parameters.   If  the  system  has  an  oscillatory 
response,  the  parameters  most  commonly  used  are  C  and  co  , 
the  damping  ratio  and  undamped  natural  frequency  respective- 
ly.  The  use  of  these  two  parameters  is  quite  natural  since 
the  coefficients  of  the  system  quadratic  equation  are  simple 

functions  of  C  and  co  . 
n 

As  seen  above,  dynamic  specifications  will  not  always 
be  given  in  terms  of  K.   and  o)  ,  but  translation  to  this  form 
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is  possible.  As  an  example,  Instrument  Engineering,  Vol.  II  -^ 
provides  graphs  of  calculated  second  order  response  in  both 

the  time  and  frequency  domains  which  may  be  used  for  this 

1^ 
purpose.  Chu   gives  several  rules  for  translating  time 

9 
domain  specifications  into  C  and  co  and  Thaler  gives  others. 

Thaler  also  includes  simplified  graphs  giving  relationships 

between  frequency  domain  specifications  and  C  and  (x>     as 

well  as  rules  of  thumb  relating  phase  and  gain  margin  to 

height  of  resonant  peak. 

Translation  of  specifications  into  C  and  co  or  relatable 
parameters  may  be  other  than  precise.  Analytical  relation- 
ships may  not  exist,  as  in  the  case  of  phase  and  gain  margins. 
Specifications  may  put  varying  demands  on  the  physical  com- 
ponents of  the  system  such  that  time  dependent  coefficients 
are  required.   They  may  prescribe  a  performance  unrealizable 
by  an  order  lower  than  third.   In  cases  such  as  these, 
reversion  to  trial-and-error  design  techniques  may  be  dic- 
tated.  If  the  response  of  a  higher  order  system  can  be  made 
to  closely  approximate  a  second  order  system,  then  adequate 
correlation  between  response  and  the  C,  o)  type  specifications 
is  possible.   It  is  to  the  problem  of  assuring  dominance  of 
a  given  quadratic  factor  that  we  now  turn  our  attention. 

A  quadratic  factor  is  said  to  have  "dominance"  if  the 
response  to  the  higher  order  system  closely  approximates  that 
of  the  specified  quadratic.  A  high  "degree  of  dominance" 
means  a  close  approximation,  i.e.  small  deviation  of  actual 
from  second  order  response.   The  use  for  which  the  system 
is  intended  and  the  strictness  of  the  specifications  will 
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indicate  the  degree  of  dominance  required. 

In  an  oscillatory  second  order  system,  the  correlation 
between  coefficients  of  the  characteristic  equation  and 


the  C  and  co  specifications  is 
n 

s^  +  2Ca)„s  +  0)   =0 
n     n 


(^.16) 


If  our  system  has  a  cubic  characteristic  equation  with 
one  pair  of  oscillatory  roots  it  will  appear  as  follows  on 
the  s -plane: 


o- 


<r-  ^^, 


1  n 


0) 


cr-> 


Figure  ^-1 

where  a,,  the  dominance  factor,  is  defined  as  shown  in 
figure  ^-1.  The  model  characteristic  equation  in  terms  of 
a-,,  C  and  0)  is: 

(s+a^a)^)(s^+2Ca3^s+co^^)  =  0 

s-^+(2C*-a3^)a)^s^<-(2Ca^+l)co^^s+a^a)^^  =  0   (h.ll) 

The  degree  to  which  the  response  of  this  system  approximates 
a  second  order  response  depends  on  the  size  of  a-j^,  A  detailed 
examination  of  the  effect  of  a,  on  the  deviation  of  a  third 
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order  response  from  that  of  a  second  order  is  shown  in 
Appendix  A. 

A  system  with  a  fourth  order  characteristic  equation 
might  take  either  of  the  following  s-plane  configurations: 
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Figure  '4-2a 
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Figure  ^-2b 


The  model  characteristic  equation  in  the  case  shown  in 
figure  ^-2a  is: 

(s+a,a)^)(s+a^a)^)(s^+2Cco„s+GO„^)  =  0 
in     d  n  n   n 

s^+(2C+P-L)a)j^s^+(l+2Cp-L+P2^^n^s^'*"^Pl"^2^p2^^n^^"*'^2^n  ^  ° 


For  the   case   in  figure  M--2b  the  model   characteristic 
equation  is: 

(s^+2aTCo^„s+aT^co  ^)(s^+2CtC0   s+co  ^)    =    0 
l^nln  Inn 

s^+(2C;j^+2a^C2)^n^^"^^^'*"^''l^l^2"^''l^^^n^^^"^^^''l^2"*"^^l''l^^'*^n^^ 
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In  choosing  the  a's  for  the  case  in  figure  ^-2a  a 
reference  to  the  curves  of  Appendix  A  is  helpful  but  not 
exact.  Each  a  must  be  sufficiently  large  so  that  its 
effect  alone  would  be  within  the  tolerances  allowed.  A 
conservative  choice  of  the  individual  a's  is  necessary  to 
insure  against  the  possibility  of  accumulating  error  of  the 
same  sign.   Determination  of  the  effect  of  more  than  one 
real  root  is  a  fruitful  area  for  further  work. 

In  the  case  where  there  are  two  pairs  of  oscillatory 
roots,  there  is  a  slightly  different  problem.   Here  it  is 
not  only  the  decay  of  the  non-dominant  pair,  but  the  size 
of  the  coefficient  of  their  sinusoidal  contribution  to  the 
time  response,  which  determines  their  effect.   In  the 
extreme  case,  where  they  are  on  the  imaginary  axis,  this 
coefficient  is  approximately  l/a  times  that  of  the  domi- 
nant pair. 

The  extreme  condition  described  above  also  serves  to 
place  a  lower  limit  on  the  size  of  the  coefficients  of  the 
characteristic  equation  in  order  to  have  stability.   If 
Cp  =  0  the  minimum  value  of  the  coefficients  of  the  fourth 
order  characteristic  equation  are  defined  as  follows: 

s^+(2CT)a)„s^+(l+aT^)co  ^s^-^{2t:.a^^ujs+a^^ai^^   =  0   (^.20) 
In       In       lln    in 

Similar  examination  of  higher  order  characteristic 
equations  shows  the  relationship  between  a's,  C  and  co  , 
and  their  coefficients.  Appendix  B  lists  the  various  forms 
of  these  coefficients,  including  minimum  stability  values, 
for  model  equations  up  to  seventh  order. 
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As  a  result  of  the  considerations  of  static  and  dynamic 
specifications  we  can  now  write  model  characteristic  equations 
for  a  system  which  will  exhibit  a  specified  response.  Con- 
siderations of  chapters  II  and  III  have  given  us  an  actual 
characteristic  equation  which  must  be  forced  into  this 
specified  configuration  by  compensation. 
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CHAPTER  5 

\ 
EFFECT  OF  COMPENSATION  ON  THE  COEFFICIENTS 

OF  Tiff;  CHARACTERISTIC  EQUATION 

5.1   Introduction*   In  Chapter  3  we  formed  the  characteristic 
equation  which  represents  the  components  in  our  servomechanism. 
The  static  and  dynamic  specifications  were  related  to  coeffi- 
cients of  model  characteristic  equations,  and  these  are  tabu- 
lated in  Appendix  B.   The  third  phase  of  our  compensation 
procedure  consists  of  determining  the  effect  of  compensation 
on  the  coefficients  of  the  characteristic  equation.   This  effect 
can  be  considered  under  two  headings: 

a)  The  Effect  of  Compensator  Location 

b)  The  Effect  of  Compensator  Type. 

The  purpose  of  this  chapter  is  to  show  how  to  quickly  deter- 
mine each  of  these  effects.   The  relationship  between  them 
and  the  coefficients  will  be  examined.  A  notation  will  be 
introduced  which  permits  rapid  evaluation  of  the  effect  of 
all  possible  compensators,  without  repeated  expansions  of  the 
characteristic  determinant. 

The  separation  of  the  two  effects  is  possible  largely 
because  of  the  unique  manner  in  which  the  addition  of  a  com- 
pensator alters  the  characteristic  determinant.   Before  con- 
sidering the  two  effects  separately  we  will  examine  the  effect 
of  compensation  on  the  determinant  of  a  general,  four-node 
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uncompensated  servomechanism. 

5.2  The  Effect  of  Compensation  on  the  Characteristic  Deter- 
minant.  In  studying  the  effect  of  a  compensator  on  the  charac' 
teristic  determinant,  we  will  use  the  four  node  system  in 
figure  5.1. 


->©- 
■^ 


aA 


-^^K^ 
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'cC 


-^- 


'dD 


•Da 

Figure  ^-1 
The  characteristic  determinant  for  this  system  is: 

0 


A 


aA 


bB 


cC 


^dD^a 


The  characteristic  determinant  is  single  valued,  and  so  its 
value  is  independent  of  the  method  used  in  expanding  it. 
For  exEimple:   When  expanded  by  cof actors  of  the  fourth  column 
5, 


A=   (-D'C-G^Sa^ 


-G 


aA     "    ^ 

°   -Sb   1 

0      0    -G 


cC 


+(-1)^(1) 


-^aA^ 


SbI 
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A  =  GclD°Da°aA°bB°cC*l 

A  similar  expansion  by  cofactors  of  the  third  row  yields: 


A  =  (-l)^(-GbB) 


1 
■^aA 
0    -G 


0  ^dD^Da 


cC 


+  (-iru) 


aA 


1   0 
0   1 


A  -   G,  r,G,r^G^^G^pG 


bB  dD  Da  cC  aA 


+  1 


Note  that  both  expansions  give  the  same  result,  and  indeed, 
any  expansion  would  give  this  result. 

If  we  were  to  add  a  compensator  to  any  path  previously 
open,  we  would  put  a  transfer  function  in  a  determinant  posi- 
tion previously  zero.   As  an  example,  consider  closing  the 
path  from  the  output  of  block  G,  g  to  node  "a"  with  an  element 
whose  transfer  function  is  Gg  .   The  characteristic  deter- 
minant of  this  modified  system  is: 


1 

^bB^Ba 

0 

^dD%a 

A„  = 

-^aA 

1 

0 

0 

0 

-^bB 

1 

0 

0 

0 

-G.p 

1 

Expanding  by  cofactors  of  the  first  row: 
10   0 


A„,  =  (-1)  (1) 
* — ^m 


-G 


bB   ^ 
0   -G 


cC 


^(-l)^(G^BSa^ 


-G 


aA 
0 


0 

1 

-G 


cC 


^0 


An 


(^dD^a^ 


1  ^  ^bB^a^aA 


aA 
0 


-G 


bB 


0 

1 


^dDSa^aA^bB^cC 


This  is  merely  the  old  A  with  the  term  G^gGg^^G^^  added. 

A  similar  effect  is  observed  is  an  existing  path  is 
paralleled.   Suppose  an  element  whose  performance  function 
is  G'   was  put  in  parallerl  with  G  «  .   This  would  make  the 
transfer  function  between  nodes  "a"  and  "b"  equal  (G  .  +  G'  ) 
Since  no  new  position  in  the  determinant  was  filled,  the 
expansion  remains: 


A,  =  (G^A 


G^A^^bB^cC^dD^a  ^  1 


An  "  ^aA^bB^cC^dD^a  ""  ^  ""  ^aA^bB^cC^dD^a 

This  again  is  the  old  /\   with  a  new  term  added.  We  can  gene- 
ralize the  above  observations  as  follows: 

ZX„  =  A  -  G^jAij  (5.1) 

where   G 


ij 


is  the  element  added  to  the  i   row,  j!   column  of 


the  determinant  as  a  result  of  adding  the  compensator.  /\ 


ij 


is  the  signed  minor  (cofactor)  of  this  determinant  position, 


It  is  the  effect  of  G.  ./X-t  ^  on  the  coefficients  of  the  charac- 
1 J   -'-J 

teristic  equation  that  we  must  determine. 

All  the  terms  in  equation  (5*1)  contain  ratios  of 
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polynomials  in  "s".   These  polynoiriials  come  from  the  products 
of  terms  G  jj,  G„  ,  etc.  in  the  characteristic  determinant. 
These  transfer  functions  represent  components  of  the  servo- 
mechanism.   By  their  physical  nature,  components  normally 
contribute  a  time  lag  to  the  system.   This  lag  appears  in 
the  transfer  function  as  an  (s+a)  factor  in  the  denominator. 
The  expansions  of  /\   and  A  consist  of  unity  plus  a  series 
of  products  of  such  transfer  functions.   Because  of  the  pre- 
valence of  lagging  components  these  products  will  have  nume- 
rators of  order  less  than  (or  at  best  equal  to)  their  denom- 
inators. When  the  terms  comprising  /\   are  placed  over  a  common 
denominator,  the  numerator  so  formed  is  the  usual  form  of 
the  characteristic  equation.   The  common  denominator  will  be 
the  product  of  all  the  non-repeated  (s+a)  factors,  and  the 
highest  power  of  s  in  any  of  the  individual  product's  denom- 
inators.  The  numerator  contains  unity  times  this  common 
denominator.   It  is  this  contribution  to  the  numerator  that 
contains  the  highest  power  of  s,  and  thus  sets  the  order  of 
the  characteristic  equation. 

With  these  thoughts  in  mind,  consider  again  equation  (5»1)! 
An  =  A  -  G^Aj  '5.1) 

Modified  Char.Egn.   _  Char .Egn.  .    n      /\  fc  o\ 

Modified  Char.  Eqn.  =  [£l£:ilm]  Char.  Eqn,  +  G.  -A^  ^CD  .  ^ 

(5.3) 
From  this  we  see  two  types  of  changes  to  the  coefficients  of 
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the  characteristic  equation: 

a)  If  rr^  is  not  unity,  i.e.,  if  the  compensator 
added  new  (s+a)  factors  or  a  higher  power  of  s  to  the  common 
denominator,  the  old  characteristic  equation  will  be  multi- 
plied by  the  product  of  these  added  terms.   This  will  raise 
the  order  of  the  characteristic  equation. 

b)  The  coefficients  of  s  in  ^i^Ai-jCDv^   will  be 
added  to  the  coefficients  of  s   in  fcDz^^ml  characteristic 
equation. 

Had  the  purpose  of  this  chapter  been  merely  to  show  the 
effect  of  a  compensator  on  the  coefficients  of  the  charac- 
teristic equation,  it  would  be  fulfilled.   One  would  only 
have  to  apply  equation  (5»3)  for  each  compensator-path 
combination.   For  an  "n"  node  system  and  "m"  possible  com- 
pensators, this  requires  (n  x  n  x  m)  applications  of  equation 
(5.3) •   This  is  a  process  that  requires  many,  many  hours  of 
algebra.   It  usually  just  isn't  done.   The  designer  relies 
on  his  previous  experience  to  help  him  pick  a  usable  path  and 
is  likely  to  overlook  one  which  would  have  done  the  job 
easier.   The  remainder  of  this  chapter  will  develop  a  tech- 
nique for  examining  all  possible  paths  in  a  reasonable  time, 
and  for  displaying  this  information  in  a  useful  manner. 
5.3  The  Effect  of  Compensator  Location.  We  have  seen 
(equation  5*1)  that  adding  a  compensator  has  the  effect  of 
adding  a  term,  G . ./\. . ,  to  the  characteristic  determinant. 
Factoring  the  performance  function  of  the  compensator  from 

the  term  G. ./\,  .   will  leave  the  effect  of  closing  the  path 
1 J   1 J 
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in  which  the  compensator  was  inserted.   It  is  this  effect 
that  we  call  the  "Effect  of  Compensator  Location".   '-'L~'^i-|« 

It  is  the  term  which  would  appear  as  characteristic  deter- 
minant element  i-j,  if  the  path  represented  by  that  element 
were  merely  closed. 

Rewriting  equation  (5.3)  for  this  condition: 

Modified  Char.  Eqn.  =  §§^  Char.  Eqn.  +  ^c^j/X^^^^^ 

(5A) 

Since  there  is  no  denominator  in  the  performance  function 
Gp  =  1,  there  are  no  new  factors  added  to  the  common  denomi- 
nator by  the  inclusion  of  G^G^/\^..   CD/v^  =  CD^  • 

Modified  Char.  Eqn.  =  Char.  Eqn.  +  ^L^ii^^A   ^^'5) 

Modified  Char.  Eqn.  -  Char.  Eqn.  =  G^y^.  .CD^   (5.6) 

(Modified  Char.  Eqn.  -  Char.  Eqn.)  will  be  a  polynomial  in  s 
containing  only  those  powers  of  s  whose  coefficients  were 
changed  by  closing  the  path.   If  we  could  quickly  determine 
the  powers  of  s  in  Gj^A .  -CDa  ,  for  all  possible  paths,  we 
would  know  what  coefficients  of  the  characteristic  equation 
were  altered  by  closing  those  paths. 

To  accomplish  this  we  will  introduce  a  new  notation. 
This  notation  will  describe  the  powers  of  s  contained  in  the 
polynomials  which  comprise  the  terms  in  equation  (5*6) . 
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Since  these  polynomials  in  s  arise  from  the  multiplication 
or  division  of  factors  (s+a)  or  "s",  they  will  contain  all 
the  successive  powers  of  s  between  the  highest  and  lowest 
power  present.   The  symbol  so  used  will  be  called  an  "Order 
Number",  and  will  consist  of  two  numbers,  one  Arabic  and 


one  Roman,  i.e.  (A,R).   The  highest  power  of  s  in  the  poly- 
nomial  is  s   ,  and  the  lowest  power  of  s  is  s  .  All  powers 
of  s  between  these  are  present.   For  example: 

If  (Order  of  P(s))  =  (3,1),  then  P(s)  =  a^^s  +a^s-'+a2S  +a^i 
The  algebra  of  these  Order  Numbers  is  much  like  that  of 
exponents.  When  two  polynomials  are  multiplied,  their  Order 
Numbers  add  to  give  the  Order  Number  of  the  product.   If 
one  polynomial  in  s  is  divided  by  another,  the  difference 
in  their  Order  Numbers  is  the  Order  Number  of  the  quotient. 
Order  Numbers  add  and  subtract  as  follows: 
(A-L,R^)  +  (A2,R2)  =  (A-lIA^,  R-LtR2) 

As  an  example  of  this  use  of  Order  Numbers,  consider: 
(s+a)(s  +bs)  =  s"^  +  (a+b)s  +abs 
Order  Number  of  (s+a)  =  (1,0) 
Order  Number  of  (s  +bs)  =  (1,1) 

(1,0)  +  (1,1)  =  (2,1) 

^         2 
Order  Number  of  s"^  +  (a+b)s  +abs  is  (2,1) 

Since  the  terms  of  equation  (5o6)  are  polynomials  in  s  of 

the  type  representable  by  Order  Numbers  we  can  rewrite  this 

equation  as  follows: 

(Coef.  Indicator)  =  (Order  of  G^^)  +  (Order  of  A-i^)  "*• 

(Order  of  CD  a  )     (5.7) 
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Before  examining  the  use  of  equation  (5*7),    it  is  well 
to  become  more  familiar  with  these  Order  Numbers,  and  to  see 
how  they  are  found  for  the  terms  of  the  equation. 

The  Arabic  number  equals  the  highest  power  of  s  arising 

from  the  combination  of  the  terms  of  the  form  (s+a). 

Example  1:   In  (s+a) (s+b) (s+c) ,  the  Arabic  number  is  3» 

Example  2:   In  Xs±aJjLs±b2>  ^^e  Arabic  number  is  1. 
(s+c) 

The  Roman  number  equals  the  net  power  of  s  arising  from 

terms  of  the  form  s  . 

Example  3:   In  s  ,  the  Roman  number  is  III. 

Example  *+:   In  s^,  the  Roman  number  is  (m-n). 
s 

As  examples  of  polynomials  in  which  both  Arabic  and  Roman 

numbers  occur: 

Example  5:   The  order  of  (s+a)   is  (0,1) 

sTs+b) 

Example  6:   The  order  of  s^(s+a) (s+b) (s+c)  is  (3,11). 

Example  7:   The  order  of  0  is  (-co, -co).   To  simplify  the 

notation,  this  will  be  designated  (-,-).   Note  that  (A,R)+ 

(-,-)  =  (-,-).   Any  sum  containing  (-,-)  has  no  value.   It 

is  as  if  you  multiplied  by  zero. 

The  Order  Numbers  comprising  the  left  side  of  equation 
i5»7)   are  obtained  as  follows: 

(Order  of  Q^):   In  determinant  positions  corresponding 
to  feedback  paths,  G^  is  the  transfer  function  of  the  block 
immediately  preceding  the  feedback  pickoff .   It  would  be  of 
the  form  of  Example  5  above.   In  determinant  positions  corres- 
ponding to  feed  forward,  or  parallel  to  a  direct  path,  it 


he 


equals  unity.   The  Order  Number  of  a  constant  is  (0,0), 

(Order  of  A.  ..);   This  could  be  found  by  expanding  the 
minor,  A . ...  and  then  applying  the  rules  exemplified  above. 
A  simpler  method  consists  of  replacing  each  existing  term 
in  the  .^characteristic  determinant  with  its  Order  Number, 
The  minors  are  then  expanded  by  the  usual  determinant  rules 
with  the  following  exceptions: 

1)  The  Order  Numbers  are  added  where  the  terms  they 
represent  would  have  been  multiplied, 

2)  Arabic  and  Roman  numbers  are  added  separately. 

3)  The  largest  Arabic  and  the  largest  Roman  numbers 
resulting  from  these  additions  comprise  the  Order  Number  of 
the  minor. 

(Order  of  CDa  ):   The  common  denominator  of  /\  will  have 
already  been  found  in  the  determination  of  the  characteristic 
equation  (Chap.  3).  Apply  the  rules  above  to  this  denomi- 
nator as  in  Example  6. 

With  the  use  of  the  Order  Number  technique,  a  matrix, 
showing  the  results  of  equation  {5»7)   for  each  path,  can 
be  formed.  This  "Location  Matrix"  is  constructed  as  follows: 

a)  Replace  the  existing  transfer  functions  in  the  Char- 
acteristic determinant  by  their  Order  numbers. 

b)  Evaluate  the  (Order  of  Aj^..)  for  each  i-J  position 
usable  for  compensation.   (Normally  all  positions) 

c)  To  this  (Order  of  A^^)  acid  (Order  of  Qj)   for  the 
position  concerned. 

d)  Add  (Order  of  G^^.)   to  (Order  of  CDa  )  and  display 
the  resulting  Order  Numbers  in  a  matrix.  Each  Order  Number 
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will  occupy  a  position  in  the  matrix  corresponding  to  the 
position  of  its  path  in  the  Characteristic  Determinant, 

This  matrix  will  indicate  which  coefficients  can  be 
changed  by  closing  a  given  path. 

To  clarify  the  ideas  presented  in  this  section,  the 
location  matrix  for  a  sample  problem  will  be  formed.   Con- 
sider the  system  of  Figure  5-2. 


Input^ 


K. 


->®- 


(s+d) 


Is+cl-^KcH   K 


ilirly 


-^-^  \ 


s(s+a 


Output  ^ 


Figure  5-2 
This  system  has  the  following  characteristic  determinant: 


A 


-K 


sTs+aT 
0 


-(s+c 
s+d 


1 

-^2 
Ts^bT 


Expanding  this  in  terms  of  elements  of  the  fourth  column: 

K-LK2K^^(s+c) 
A-     ^  ^   s(s+a)(s+b)(s+d) 


A  _  s(s+a)(s+b)(s+d)+K-^K2K^(s+c) 
s(s+a)(s+b)(s+d) 
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Expanding  the  numerator  and  introducing  standard  symbols  for 

coefficients: 

^32 
/-.v.   -m        s  +a-,s  +a^s  +a-,  s+a 
Char.Eqn,   _     ^  2    1   o 

CD  A  s(s+a)(s+b)(s+d) 

Applying  the  rules  for  formation  of  a  Location  Matrix: 

(0,0)    (-,-)     (-,-)     (-1,-1) 

(0,0)    (0,0)     (-,-)       (-,-) 

(-,-)    (0,0)     (0,0)       (-,-) 

(-,-)    (-,-)      (-1,0)      (0,0) 


Step  a: 

A 


step  b:   All  ■  ^°'°^  ^12  "  ^°'°^  ^13  '  ^°'°^  ^l^f  "  ("IjO^ 

A21    -    (-2,-1)  A22   -    (0»0)  A23   -   (0,0)  A2I,   -   (-1,0) 
A31    -    (-2,-1)   A32    -    (-2,-1)   A33    -    (0,0)   A3I,   -   (-1,0) 


Ai,i  -  (-1,-1)  ZV2  -  (-1,-1)  Ai,3  -  (-1,-1)  /\i 


(0,0) 


The   Order  numbers   of   the   G^-'s   are   as   follows: 

C  C*  C  P 

Hi  -  (0,0)     %2  -  (0,0)     ^3  -  (-1,0)     hh  -  (-1,-1) 

err  c 

hi  -  (0,0)     ^2  -  (0,0)     ^3  -  (-1,0)     hh  -  '(-1,-1) 
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(0,0)     ^2   -   (0,0)      h3   -   (-1,0)      h^  -   (-1,-1) 


c  r 

\l    -    (0,0)      \2    -    (0,0) 


c  c 

\3  -  (  0,0)     ^kh  -  (-1,-1) 


step  c:   (Order  of 

OiAij)  -  po 

sition  in  the 

array  denotes 

subscripts. 

(0,0) 

(0,0) 

(-1,0) 

(-2,-1) 

(-2,-1) 

(0,0) 

(-1,0) 

(-2,-1) 

(-2,-1) 

(-2,-1) 

(-1,0) 

(-2,-1) 

(-1,-1) 

(-1,-1) 

(-1,-1) 

(-1,-1) 

(Order   of  CD  a  )   equals    (3,1) 
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Step  d: 


Location  Matrix 


(3,1) 

(3,1) 

(2,1) 

(1,0) 

(1,0) 

(3,1) 

(2,1) 

(1,0) 

(1,0) 

(1,0) 

(2,1) 

(1,0) 

(2,0) 

(2,0) 

(2,0) 

(2,0) 

Step  b,  finding  the  Order  of  each  minor,  merely  entails 
scanning  the  auxiliary  determinant  of  step  a  to  find  the 
largest  Arabic  and  the  largest  Roman  sums  present.   The 
presence  of  numerous  (-,-)  terms  in  an  uf^compensated  system 
makes  this  a  few  minutes  work.   This  is  in  comparison  to 
the  hours  of  algebra  required  to  evaluate  0.  ./\.  .CDa  for 

every  possible  compensator. 

The  Location  Matrix  is  the  sum  of  the  terms  of  Equation 
{5 '7)   for  all  possible  compensation  paths.   The  information 
in  the  Location  Matrix,  modified  by  the  effect  of  compensator 
type,  will  indicate  changeable  coefficients  of  the  Charac- 
teristic Equation. 

5.3 •   The  Effect  of  Compensator  Type 
5.3-1   Introduction 

Examination  of  the  type  of  compensator  yields  informa- 
tion on: 

a.  Which  coefficients  of  the  Characteristic  Equation 
will  be  changed,  and 

b.  the  extent  of  that  change. 

It  is  the  Order  Number  of  the  compensator  in  conjunction 
with  the  Location  Matrix  that  indicates  which  coefficients 
will  be  changed.   The  performance  function  indicates  the 
magnitude  of  the  change. 
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Table  V-1  lists  the  compensator  types  to  be  considered, 
with  their  Order  Number  and  Performance  Function. 


Table  V-1  Elementary  Compensators 

TYPE 

ORDER  NUMBER 

PERFORMANCE 
FUNCTION 

Sensitivity 

(0,0) 

K 

Differentiator 

(0,1) 

s 

First  Order  Lead  Filter 

(1,0)' 

(s+a) 

Second  Order  Lead  Filter 

(2,0) 

(s^+bs+c) 

Integrator 

(0,-1) 

1/s 

First  Order  Lag  Filter 

(-1,0) 

l/(s+a) 

Second  Order  Lag  Filter 

(-2,0) 

l/(s^+bs+c) 

5.3.2  Determination  of  Which  Coefficients  are  Changed 

The  Order  Number  of  the  compensator,  when  added  to  the 
element  of  the  Location  Matrix  corresponding  to  its  path, 
gives  a  new  Coefficient  Indicator.  As  previously  defined, 
the  (A,R)  pair  resulting  from  this  addition  indicates  which 
coefficients  of  the  Characteristic  Equation  are  modified; 

There  are  three  situations  in  which  addition  of  a  com- 
pensator will  raise  the  order  of  the  Characteristic  Equation. 
If  the  order  is  raised,  there  will  be  two  changes  to  the 
coefficients.   One  change  will  be  to  all  coefficients  due 
to  merely  raising  the  order.   The  change  dictated  by  the 
Coefficient  Indicator  will  be  an  additional  change  to  the 
coefficients  of  the  higher  order  equation. 

Before  examining  the  quantitative  effects  of  compensation 
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on  the  coefficients,  the  three  cases  where  the  order  of  the 
Characteristic  Equation  is  raised  will  be  considered: 

Case  I  -  This  situation  arises  when  integrator  com- 
pensation is  used.  If  the  addition  of  a  com- 
pensator results  in  one  of  the  terms  of 

l   +  EL+Z7tL^  +  E7rL+ 

having  more  integrations  than  the  overall 
uncompensated  closed  loop  performance  func- 
tion, the  order  of  the  Characteristic  Equation 
will  be  raised.   In  the  above  expression,  L 
signifies  loop  performance  function 


TT 

k 


signifies  products  of  loop  performance 
functions  taken  k  at  a  time.   That  the  order 
of  the  Characteristic  Equation  is  raised  in 
this  situation  is  apparent  if  the  Character- 
istic Equation  is  derived  from  the  system  per- 
formance function  obtained  by  application  of 
Mason's  Rule. 

A  simple  example  of  this  effect  may  be 
seen  by  considering  the  following  system: 
Block  A  Block  B 

^ — 


Kn 


K, 


52 


Where  A  = 

1 

s 

A  = 

1  -» 

1 

-  ^1^2 
s 

Char.  Ean. 

_  s_ 

+  ^1^2 
s 

Char.  Eqn.  =  s  +  V2   ° 


For  this  system  the  Location  Matrix  is: 

(0,1)      (0,0) 
(0,0)       (0,0)^ 

Consider  integration  feedback  from  the  output  of 
Block  B  to  Node  b: 


L.  M, 


Here,  A 


Block  A. 

^^     ] 

Block  B 

4 

^y-^ 

% 

K2 

^ 

^ 

^     ^ 

/ 

\ 

s 

1 
s 

's 

A  = 


■K.      1+  ^2 

1        "~2" 
s 


1   +  ^  -»-  ^1^2 
s^    ^ 
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Char.  Ean.  _   s^  +   KiKps  +   K2 

2 
Char.  Eqn.  =  s  +  K^K2S  +  K^  =  0 


Notice  that  the  order  of  the  Characteristic 
Equation  has  been  raised  by  this  combination 
of  compensator  and  path.   If  the  integration 
had  been  fed  back  from  the  output  of  Block  A 
to  Node  a,  a  loop  performance  function  higher 
than  first  order  would  not  have  been  generated 
and  the  order  of  the  equation  would  not  have 
been  raised.   In  this  case, 


A  = 

=   1 

s               s 

Char.   Eqn.   _ 

s_ 

+    (Ki)(l+K^) 
s 

Char.    Eqn.    = 

s  +  K^(l+K2)   =   ( 

Table  V-1  gives  an  Order  Number  of  (0,-1) 
for  this  compensator.  Adding  this  to  the 
proper  element  of  the  Location  Matrix  gives; 

Fed  from  B  to  b:  (Coef.  Ind.)  =  (0,-1) 
Fed  from  A  to  a:  (Coef.  Ind.)  =  (0,0) 
Since  the  Characteristic  Equation  contains  no 
negative  powers  of  s,  the  first  indicator  (for 
B  to  b)  appears  to  be  contradictory.   It  is 
this  very  contradiction  that  yields  the 
information  on  equation  order.   The  order  of 
the  Characteristic  Equation  must  be  raised 
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by  the  amount  that  the  "R"  of  the  Coefficient 
Indicator  is  negative.   That  this  occurred  is 
borne  out  by  comparison  of  the  Characteristic 
Equations  in  the  above  example  before  and  after 
this  compensation.   In  the  case  of  feedback  from 
A  to  a,  the  "R"  of  the  Coefficient  Indicator 
accurately  reflects  that  there  is  no  change  in 
order. 

Rule  I  -  The  order  of  the  Characteristic  Equation  is 
raised  by  the  amount  that  the  Roman  Numeral  of  the 
Coefficient  Indicator  is  negative. 

Case  II  -  Another  order  raising  situation  arises  when 
any  of  the  terms  of 

1  +  Z  L  +  ZTTL  +  ZTT  L  + 

r^   T   2     ^        T   3 

when  placed  over  a  lowest  common  denominator,  has 
a  numerator  of  higher  order  than  the  uncompensated 
equation.   Again,  this  is  apparent  if  the  Charac- 
teristic Equation  is  derived  from  the  system  per- 
formance function  obtained  by  application  of 
Mason's  Rule. "'"■'" 


For  example,  if  the  sample  system  is  compen- 
sated as  follows: 

Block  A  Block  B 


r\ 



h> 

K2 

\ 

r 

2 
s  +bs+c 
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then,  A  = 


-K^   1+  K2(q  +bs+c) 


J 


A  =  1  -^  K1K2  +  K2(s  +bs-^c) 
s  s 

Char.  Eon.  ^  K2S^-^  (K2b'H)s  +  (K2C  +  K1K2 


CD 


A 
Char.  Eqn. 


(k^ 


(K^+c) 


Adding  the  Order  Number  of  this  compensator,  (2,0), 
to  the  proper  Location  Matrix  element  gives  a 
Coefficient  Indicator  of  (2,0).  Again  there  is  an 
apparent  contradiction  in  that  the  highest  coeffi- 
cient indicated  to  be  changed  is  not  present  in 
the  uncompensated  Characteristic  Fquation.   This 
contradiction  indicates  that  the  order  of  the  Char- 
acteristic Equation  must  be  raised. 

Putting  this  compensator  in  the  path  from  the 
output  of  Block  A  to  Node  a  yields: 
A  ^  Char.  Eqn.  ^  K^s^  -^  Kibs^  +    (Kic+l)s  +  K1K2 

Char.  Eqn.  =  K-^s^  +  K^bs^  +  (K^c+l)s  +  K-^K2  =  0 

The  Coefficient  Indicator  in  this  case  is  (2,1), 
which  again  indicates  a  raise  in  the  order  of  the 
Characteristic  Equation,  the  coefficients  to  be 
changed  being  a^?  a2  and  a^* 
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Rule  II  -      The  order  of  the  Characteristic  Equation  is 
raised  by  the  amount 

(A^K^Coef.  Ind.  "  ^^*^^Oviev   of  CD^ 

Note:   If  an  occasion  arises  where  both  Rules  I 
and  II  apply,  first  raise  the  Order  of  CD  a  by 
the  amount  indicated  in  Rule  I,  put  "R"  of  the 
Coefficient  Indicator  to  zero,  then,  if  necessary, 
apply  Rule  II. 

Case  III  -      A  lag  filter  will  always  introduce  its  denom- 
inator in  CD.  .   This  will  raise  the  order  of  the 
common  denominator  by  an  amount  equal  to  the  order 
of  the  filter.   Since  in  practical  servos,  the 
order  of  the  Characteristic  Equation  equals  the 
order  of  the  common  denominator  of  /\.,  it  too 
will  be  raised  by  this  amount. 

Rule  III  -  The  order  of  the  Characteristic  Equation  is 
raised  by  the  amount  the  Arabic  Numeral  of  (Order 
of  the  Compensator)  is  negative. 

5. 3*3  Determination  of  Quantitative  Changes  to  Coefficients 

The  Modified  Characteristic  Fquation,  Equation  (5«3), 
consists  of  two  terms: 

Mod.  Char.  Eqn.  =  '^^Am  (Char,  Eqn.)  +  G.  -A^  ^CDa^ 
CDa  ■^*'  ■^♦^  ^-^^ 

This  indicates  that  there  are  two  contributions  to  the  change 
in  the  coefficients.   The  first  term's  effect  is  predictable 
in  the  case  of  simple  compensators.   It  is  the  amount  all 
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coefficients  are  changed  when  the  order  of  the  Characteristic 
Equation  is  raised  by  applying  Rules  I  and  III  above. 


In  Case  I,   CDy^^   g|R|   where  the  "R" 


*^^A        is  the  Roman 

Numeral  of  the 

Coefficient 

Indicator. 

In  Case  III,  ^^Am  ^  (s+a)  or  (s^+bs+c) 
CD/^  - 

*— ^       depending  on  the 

order  of  the  fil- 
ter used. 

In  these  cases,  the  initial  effect  on  coefficients  can  be 
readily  calculated. 

The  change  to  coefficients  a.^p>,  ^.^.^d,  etc.,  as  indi- 
cated by  the  Coefficient  Indicator  can  only  be  found  by  eval- 
uating '^i -iZXj^^CD  A   for  the  path  concerned.  Evaluating  the 

change  on  coefficients  caused  by  raising  the  order  of  the 
Characteristic  Equation  in  accordance  with  Rule  II  can  only 
be  determined  by  evaluating  G.  ./\.  .CD  a  • 
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CHAPTFR  6 

SELECTION  OF  COMPENSATOR 
1.   Introduction 

The  problem  to  be  considered  in  this  chapter  is  that  of 
choosing  the  combination  of  compensator  and  location  that  can 
most  easily  change  the  characteristic  equation  of  the  system 
to  one  established  by  the  specifications.   This  requires  a 
comparison  of  the  coefficients  of  the  uncompensated  charac- 
teristic equation  to  those  of  a  minimum  model  equation.   The 
actual  coefficients  which  exceed  a  minimum  allowable  value 
will  be  used  to  select  a  working  model  equation.   The  work- 
ing model  equation  permits  calculation  of  a  compensated  char- 
acteristic equation  consistent  with  the  specifications. 

There  are  four  equations  of  interest  in  this  chapter. 
The  uncompensated  characteristic  equation  describes  the 
servomechanism  before  compensation.   Its  coefficients  are 
set  by  the  basic  hardware.   It  is  these  coefficients  which 
must  be  changed.   The  minimum  model  equation  is  the  charac- 
teristic equation  of  a  system  with  a  transient  response  which 
has: 

a.  A  specified  dominant  damped  oscillatory  mode, 

b.  A  minimum  number  of  first  order  modes,  (none  or  one), 

c.  Non-dominant  undamped  oscillatory  modes. 

No  first  order  modes  will  exist  for  even  ordered  systems.  No 
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non-dominant  oscillatory  mode  will  exist  for  a  third  order 
system.   Its  coefficients  are  those  of  the  characteristic 
equation  of  a  system  verging  on  instability.   The  working 
model  equation  is  the  characteristic  equation  of  a  system 
with  a  transient  response  which  has  a  specified  dominant 
damped  oscillatory  mode.   Its  non-dominant  modes  will  be 
damped  oscillatory  and/or  first  order.   Its  coefficients  are 
those  of  a  characteristic  equation  of  a  stable  system.   The 
compensated  characteristic  equation  describes  the  servo- 
mechanism  after  compensation.  Depending  on  its  stage  of 
development,  its  coefficients  may  be  purely  numerical  or 
may  contain  unspecified  parameters.   Compensation  which 
changes  the  coefficients  of  the  uncompensated  equations  to 
those  of  the  compensated  equation  will  result  in  a  system 
with  the  specified  response. 

2.   Choosing  the  Minimum  Model  Equation 

There  are  three  factors  which  determine  the  choice  of 
a  minimum  model  equation.   They  are: 

1.  It  must  be  of  the  same  order  as  the  uncompensated 
characteristic  equation. 

2.  The  quadratic  made  up  of  the  specified  C  and  co 
must  be  a  factor  of  the  model. 

3.  The  dominance  factors,  a's,  must  be  sufficiently 
large  to  insure  that  the  response  of  the  model  has  the 
desired  correspondence  with  the  response  of  the  quadratic 
The  minimum  model  equations  in  Appendix  B  are  so  con- 
structed that  substitution  of  the  desired  C,  (o  and  a's 
will  give  them  the  properties  listed  above. 


60 


3,  Determination  of  the  Compen-sated  Characteristic  Equation 
The  compensated  characteristic  equation  must  have  coef- 
ficients at  least  as  large  as  those  of  the  minimum  model. 
If  its  coefficients  are  less  than  those  of  the  minlpum  model, 
the  system  must  be  unstable.   A  comparison  of  the  uncompen- 
sated characteristic  equation  to  the  minimum  model  equation 
will  show  which  coefficients  must  be  changed.   There  are 
two  possible  conditions  arising  from  this  comparison.  Either 
all  of  the  coefficients  are  below  the  minimum,  or  a  number 
of  them  are  already  sufficiently  large.   In  the  first  case 
substitution  of  non-dominant  damping  ratios  and/or  acceptable 
dominance  factors  into  the  working  models  of  the  proper  order, 
as  given  in  Appendix  B,  will  set  the  compensated  characteristic 
equation.   The  ratio  of  a-,  to  ag  in  the  compensated  charac- 
teristic equation  may  be  set  by  static  specifications.   This 
may  limit  the  choice  of  a's  and  non-dominant  C's. 

In  the  case  where  some  or  all  of  the  coefficients  are 
already  above  the  minimum,  several  factors  must  be  considered. 
These  ares 

1.  The  interdependence  of  the  coefficients  of  the  model 
equation. 

2.  The  position  of  the  satisfactory  coefficients. 

3.  The  effect  of  static  specifications. 

The  interdependence  of  coefficients  is  expressed  in  the 
working  model  equations  of  Appendix  B.   In  an  n   order  model 
the  n  coefficients  are  functions  of  the  specified  C  and  co  , 
and  of  (n-2)  dominance  factors  and/or  non-dominant  C's. 
This  gives  rise  to  n  equations  in  n  unknowns.   Two  of  these 
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unknowns,  C  and  co  of  the  dominant  quadratic,  are  set  by- 
specifications.  The  remaining  (n-2)  degrees  of  freedom 
may  be  constrained  by  fixing  either  coefficients  or  a's. 
If  any  (n-2)  coefficients  of  the  uncompensated  characteristic 
equation  are  above  the  minimum  values,  they  will  determine 
the  required  value  of  the  remaining  two  coefficients.   If 
fewer  than  (n-2)  coefficients  are  satisfactory,   the 
remaining  degrees  of  freedom  may  be  fixed  by  an  arbitrary 
choice  of  a's.   Required  values  for  the  remaining  coef- 
ficients can  then  be  calculated.   These  satisfactory  and 
calculated  coefficients  comprise  the  compensated  character- 
istic equation  which  will  require  minimum  change  in  coeffi- 
cients.  The  satisfactory  coefficients  should  first  be  com- 
pared to  the  model  which  contains  all  real  non-dominant  roots. 
If  they  are  not  sufficiently  large  to  fit  this  model,  then 
the  oscillatory  working  model  equation  must  be  used. 

The  second  consideration  arises  from  the  fact  that  all 
compensators  affect  successive  coefficients  in  the  character- 
istic equation.   As  developed  in  Chapter  5?  the  independent 
coefficient  changes  are  found  by  evaluating  G. ./\. .CDa  • 

This  will  be  in  the  form  of  Ks^(s+a) (s+b)  etc.  When 
expanded,  it  will  contain  successive  powers  of  s.   This  fact 
is  also  revealed  by  the  Coefficient  Indicator.   Because  of 
this  fact,  the  satisfactory  coefficients  used  to  determine 
the  compensated  characteristic  equation  should  be  chosen  so 
that  the  coefficients  remaining  to  be  changed  are  successive. 

The  static  specifications  placed  lower  limits  on  some 
coefficients  of  the  compensated  characteristic  equation. 
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These  limits  may  force  some  coefficients  to  be  larger  than 
the  dominance  factors  dictate.   The  coefficients  so  affected 
must  be  included  in  the  (n-2)  degree-of -freedom  considerations 
above. 

Comparison  of  the  desired  Coefficient  Indicator  and  the 
Location  Matrix  will  determine  the  type  of  compensator  required 
to  change  the  necessary  coefficients.   If  it  raises  the  order 
of  the  equation,  the  abpve  procedure  must  be  repeated  start- 
ing with  the  minimum  model  equation  of  the  higher  order.   The 
uncompensated  equation  to  be  compared  with  the  model  can  be 
formed  knowing  the  order  raising  effects  described  in 
Chapter  5« 

^.   Evaluating  the  Parameters  of  the  Compensator 

The  G. -A »  .  of  the  required  compensator  is  added  to  the 
old  characteristic  determinant.   This  results  in  a  charac- 
teristic equation  in  which  the  coefficients  requiring  modifi- 
cation contain  the  parameters  of  the  compensator.   Equating 
these  coefficients  to  those  of  the  compensated  characteristic 
equation  will  determine  the  value  of  the  parameters.   The 
compensator  must  have  sufficient  degrees  of  freedom  to 
match  the  independent  coefficient  changes  required. 
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CHAPTER  7 
COMPENSATION  DESIGN  PROCEDURE 

7.1  Introduction 

The  advantages  of  the  technique  described  In  the  previous 
chapters  are  three.   First,  It  Involves  no  graphic  presenta- 
tions as  required  by  other  methods.   Second,  It  does  not  depend 
on  the  Intuition  of  the  designer  for  Its  choice  of  compensation 
to  achieve  the  desired  result.   Finally,  It  requires  little 
of  the  extensive  algebra  Inherent  In  other  analytical  methods. 
Since  there  Is  no  Intuition  or  experimental  graphing  Involved, 
an  ordered  analytical  approach  to  any  compensation  problem 
Is  possible. 

This  method  divides  naturally  Into  four  phases,  which 
are  the  subjects  of  Chapters  3  through  6.   These  chapter 
divisions  were  made  primarily  In  order  to  explain  the  tech- 
nique.  They  also  follow  In  the  proper  logical  order  for 
application  to  the  solution  of  a  problem.   The  step-by-step 
procedure  In  this  chapter  thus  serves  as  a  summary  of  the 
paper  as  well  as  a  guide  for  the  compensation:  designer. 

7«2  Stei>-bv-step  Guide  to  Compensation 

Given:   I.   A  servomechanlsm  for  which  the  direct  path 

components  have  been  selected.   This  servo  may  have 
been  partially  compensated  or  may  just  have  unity 
feedback  from  output  to  Input. 
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II.   Static  and  Dynamic  Performance  Specifications 

to  be  met  by  the  system. 
Problem:   To  design  compensation  to  meet  the  specifications. 
Solution: 

Phase  I:   Forming  the  characteristic  equation. 

A.  Arrange  the  system  in  a  standard  block  diagram. 
COMMENT:   Combine  the  components  between  usable 
pickoffs  and  nodes  into  single  blocks. 

B.  Form  the  characteristic  determinant. 

C.  Write  the  characteristic  equation  by  setting 
the  expanded  characteristic  determinant  equal  to 
zero. 

COMMENT:   Note  the  Order  Number  of  CD  a  . 
Phase  II:   Relating  Specifications  to  Coefficients  of 
Model  Characteristic  Equations. 

A.  Express  static  specifications  in  terms  of 
restrictions  on  Rq   and  a-,  of  the  compensated 
characteristic  equation. 

B.  Interpret  dynamic  specifications  in  terms  of 
second  order  response  characteristics!  i.e.,  C  and 

00  o 

n 

Co   Choose  minimum  acceptable  dominance  factors 
(a's  and  non-dominant  C's)  to  insure  second  order 
dominance. 
Phase  Ills   Determining  the  Effect  of  Possible  Compen- 
sators on  Coefficients  of  the  Uncompensated  Character- 
istic Fquation§  i.e.,  Forming  the  Location  Matrix. 
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COMMENT:   Chapter  5  shows  the  details  of  this  phase. 
Phase  IV:   Choosing  the  Compensator. 

A.  Form  the  minimum  model  equation  using  a's,  C's  and 

00^  found  in  Phase  II. 
n 

COMMENT:   All  coefficients  of  the  compensated  equation 
must  be  greater  than  those  of  the  minimum  model  to 
insure  stability. 

B.  Choose  the  Working  Model  Equation. 

1.  Insure  compliance  with  dynamic  specifications. 
Put  sufficiently  large  a's  (Phase  II-C)  and  the  required 
C  and  0)  (Phase  II-B)  in  the  working  model  equations  of 
the  proper  order. 

2.  Insure  compliance  with  static  specifications. 
Compare  a^/a  ratio,  of  the  various  forms  of  the  working 
model  equation.   Choose  for  a  working  model  that  form 
which,  with  allowable  non-dominant  zetas,  (Phase  II-C) 
has  the  smallest  coefficients  that  meet  the  specifications 

3.  Remove  the  effect  of  static  specifications  from 
working  model  equation.   In  the  event  that  the  minimum 
model  equation  will  not  meet  static  specifications,  the 
basic  system  configuration  must  be  changed  by  addition 
of  a  component  which  reduces  following  error  to  zero. 
This  will  invariably  mean  the  addition  of  an  integration. 
Compensation  of  this  type  will  introduce  a  zero  in  the 
system  performance  function.   Set  this  zero  to  a  value 
which  will  remove  its  effect  on  the  dominant  mode  of  the 
response  (See  Figure  A-2),   The  design  procedure  must 

be  started  again  with  the  revised  system. 
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Form  the  compensated  Characteristic  equation- 

1.  Compare  the  uncompensated  characteristic  equa- 
tion to  the  working  model. 

2.  If  n-m-2  coefficients  of  the  uncompensated  equa- 
tion are  above  the  working  model  coefficients  and 
leave  successive  coefficients  to  be  altered,  use 
them  to  fill  some  of  the  n-2  degrees  of  freedom. 
Set  m  arbitrary  a ' s  or  non-dominant  C's.   Compute 
the  m+2  remaining  coefficients  of  the  compensated 
characteristic  equation.   If  no  coefficients  are 
above  those  of  the  working  model,  it  is  the  com- 
pensated characteristic  equation. 

Select  the  Compensator. 

1.  Form  the  desired  Coefficient  Indicator  from  the 
coefficients  in  C-2  above  which  must  be  changed. 

2,  Compare  this  Coefficient  Indicator  to  the  Loca- 
tion Matrix  to  find  the  simplest  possible  compen- 
sator. The  compensator  selected  must  introduce 
enough  additional  variable  parameters  so  that  the 
total  number  equals  the  number  of  coefficients  to 
be  changed  plus  the  number  of  zeros  of  the  system 
to  be  set  (See  COMMFNT) 

3o  Form  In  1  J^-^i .1  ,  which  is  the  system  perfor- 

^m 
mance  function  relating  the  actuating  input  at 

node  "n"  to  the  putput  of  block  G^j"   This 

expression  will  contain  the  adjustable  parameters 

of  the  compensator.   The  compensated  characteristic 
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equation  will  be  the  denominator  of  this  function. 

k-.        Equate  the  parametric  form  of  the  coefficients  of 

the  denominator  of  the  system  function  to  the  coef- 
ficients of  the  compensated  characteristic  equation, 
Also  set  the  parameters  of  any  numerator  zeros  to 
a  value  which  gives  a  satisfactory  zero  (i.e., 
sufficiently  large  a). 

COMMENT:   If  the  only  possible  compensator  raises  the 
order  of  the  characteristic  equation,  apply  the 
order  raising  rules  of  Chapter  5  to  the  uncompen- 
sated characteristic  equation.   This  order  raising 
will  put  a  zero  in  the  system  function.  Repeat 
phase  III  using  higher  order  model. 
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CHAPTER  8 

SAMPLE  PROBLEM 

8.1  Introduction 

The  problem  to  be  worked  in  this  chapter  is  a  modification 
of  an  analysis  problem  presented  by  Chestnut  and  Mayer.    The 
problem,  as  presented  in  Chapter  13  of  the  reference,  consisted 
of  determining  the  performance  of  a  position  control  system 
subjected  to  a  load  disturbance.   Compensation  had  already 
been  accomplished.   It  consisted  of  feedback  of  the  first  and 

second  derivatives  of  the  output  quantity.   The  signal  fed 

2 

back  was  modified  by  a  filter  of  the  form  Ks  (s^^a) » 

(s+b)(s+c)(s+d) 

The  specifications  were  estimated  from  the  performance  of  the 
system  as  compensated  in  the  reference. 

3.2  Statement  of  the  Problem 

Given.' X  The  system  to  be  compensated  has  the  primary  require- 
ment of  position  control.   However,  it  must  respond  properly, 
within  specified  limits,  when  subjected  to  a  designated 
disturbance  load  torque.   The  main  path  components  consist 
of  the  following? 
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Component 
Preamplifier 
Amplifier,  Amplidyne 

Motor  and  Load 


Disturbance  Load/Voltage 
Transducer 


Transfer  Function 


K,  not  specified 

1^00       850 
(s+31.5)   (s+17.6) 

^-,78 


s(s+9.1) 


K^  =  o855  volt/ft, lb, 


II  Static  Specifications 

Steady  State  Error  in  response  to  step     =  0*0 

Following  Error  in  response  to  ramp 

Maximum  output  in  response  to  load  torque 

Dynamic  Specifications: 

Peak  Overshoot  =  1<.2 

Settling  Time  to  iVq  of  steady  state 

in  response  to  step  =  o6  sec 

Problem?   Design  compensation  such  that  the  system  will 

conform  to  specifications. 

8.3  Solution  of  the  Problem 

Phase  I:   Forming  the  characteristic  equation. 


10  rad/r ad/sec 
10"^rad/ft.lbo 


il. 


Load  Torcrue 
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A 


1 

0 

0 

^.78 
s(s+9.1) 

•% 

1 

0 

0 

0 

-1500 
s+31. 

"5 

1 

0 

0 

0 

■^:§5^ 

1 

s+17.6 


A  -  1  4.      (6.09  X  10  )  Ki 

^  "  ^   s(s+31.5)(s+17.6)(s+9.1) 


A  -  s(s+31.5)(s+17.6)(s+9.1)+(6.09  x  10  )  Ki 
^  s(s+31.5)(s+17.6)(s+9.1) 

A  -  s^+58.2s^+1001s^+50^0s+(6.09  x  10^)  Ki 
^  ~       s(s+31.5)(s+17.6)(s+9.1) 

Char,  Eqn:  s^  +  58. 2s^  +  lOOls^  +  50^0s  +  (6o09  x  10^)K^  =  0 
(Order  of  CD  a  )  =  (3,1) 
Phase  II:   Relating  Specifications  to  Coefficients  of  Model 
Characteristic  Equation, 

A.   Static  Specifications 

lo  Steady  State  Position  Error  in  response  to  step  =  0, 
^o  =  l^l^l   f^dlAi^CD^] 


Limit 


r  '+.78    (1.275  X  10^)Ki   ^ 
Ls(s+9ol)  (s+31o5)(s+17.6) 

s(s+9ol)(s+31o5)(s+17o6)] 


a^  =  (6.09  X  10'')K-L 
The  value  of  a  is  independent  of  this  specification 
because  of  the  unity  feedback. 


2.   Following  Error  in  response  to  ramp  =  10   sec 
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Error^g  =  10  >  ,.  ^^   .^6 


(6.09  X  10" )K^ 

K  >  SOitO    ^  8.28  or  ^  <   lO"^ 

^       6.09  X  10  o 

3.  Response  to  Load  Torque  =  10~  rad/ftolb. 

k  >  Limit  ["G^^^Ai^I^CD^l 

10"^  >  Limit  r(o855)(^.78)(s)(s^-17.6)(s-H^lo'^)(s+9oin 
^  "  °  L (6.09  X  10^)(s)(s+9ol)K^ J 


10-^  >  (.8^^)Clfo78)(17.6)(^l.^) 
(6.09  X  10  )K-j_ 

%  >  3-^3 
B.  Dynamic  Specifications. 

From  reference  13,  Peak  Overshoot  of  1.2  corresponds 

to  (.6  >  C  >  •5).  We  will  choose  the  conservative 

C  =  .6. 

From  reference  10,  Error  is  reduced  to  within  1  /^ 
of  final  steady  state  value  in  four  time  constants,  and 


Settling  Time  =  h     =  .6  sec 


^^n 


co_ 


d     =  — T^  =  11 ol  rad/sec 

C.   Choosing  Dominance  Factors. 

From  Appendix  A,  ^'t  >  3  and  ^p  -  ^  ^^^   acceptable 
values  of  dominance  factors.   If  non-dominant  roots  are 
complex,  use  a  >  3  and  C2  >  3o 
Phase  III  J   Determining  the  Effect  of  Possible  Compensators. 
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(0,0) 

(-,-) 

(-,-) 

(-1,-1) 

Auxiliary  /\   = 

(-,-) 
(-,-) 

(0,0) 
(-1,0) 

(-,-) 

(0,0) 

(-,-) 
(-,-) 

(-,-) 

(-,-) 

(-1,0) 

(0,0) 

(0,0) 

(0,0) 

(-1,0) 

(-2,0) 

Order  Numbers 

(-3,-1) 

(-2,-1) 

(-1,-1) 

(0,0) 
(-2,-1) 

(-1,-1) 

(-1,0) 

(0,0) 

(-2,-1) 

(-2,0) 
(-1,0) 
(0,0) 

(0,0) 

(-1,0) 

(-1,0) 

(-1,1) 

(0,0) 

(-1,0) 

(-1,0) 

(-1,-1) 

Order  Numbers 

(0,0) 

(0,0) 

(-1,0) 

(-1,-1) 

of  Gj^'s 

(0,0) 

(0,0) 

(0,0) 

(-1,-1) 

(Order  of  CD  a  )  =  (3,1) 


Location  Matrix= 


(3,1) 
(0,0) 
(1,0) 
(2,0) 


(2,1) 
(2,1) 
(1,0) 
(2,0) 


(1,1) 
(1,1) 
(2,1) 
(1,0) 


(0,0) 
(0,0) 
(1,0) 
(2,0) 


Phase  IV:   Choosing  the  Compensator 


(1+0^^)0)^^  =  1,233 


a,  =  2r:-,^.^(^  ^  =   m-,080 
1     1  1  n      ' 
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%  =  °i\^  =  137,000 

s  +  13.33s^  +  1233s^  +  1^080s  +  137000  =  0 

Note  that  a^   and  a-,  of  the  uncompensated  system  are 

below  minimum  values. 

B.   a;^  for  minimum  model  =  .103 

ai  for  all  complex  non-dominant  roots  =  .126 

an  for  all  real  non-dominant  roots  =  .156 

"o 

Since  a-./a  from  the  model  equations  can  never  be 
less  than  .103,  the  condition  on  a-,/a  imposed  by  static 
specifications  must  be  relieved  by  making  the  following 
error  independent  of  this  ratio.   This  effect  can  be 
achieved  by  compensating  so  as  to  have  zero  following 
error.   This  implies  making  the  servo  a  velocity  control 
system.   However,  it  may  not  lose  its  position  control 
ability.   Paralleling  the  main  path  with  another  forward 
path  containing  an  extra  integration  will  result  in  both 
position  and  velocity  control.   This  effect  can  also  be 
obtained  by  placing  the  extra  integration  in  a  feed  for- 
ward path  from  the  first  to  second  nodes  thereby  utilizing 
succeeding  integrations  in  both  forward  paths  simultaneous- 
ly. Using  this  compensation,  the  transfer  function,  G^j 
becomes      ^iSULl^/^ 

In  addition  to  raising  the  order  of  the  characteristic 
equation,  a  zero  of  the  system  performance  function  will 
be  introduced  at  s  =  -K2/K-,  .  Set  this  zero  such  that  its 


7^ 


a  >  2,  s^y  Kp/K^  =  25.   It  is  now  necessary  to  reconmence 

the  solution. 

Phase  I:   Forming  the  characteristic  equation, 

A  Load  Torque 


-^ 


s 


855 


^*- 


1500 
s+31.5 


-O- 


850 


s+17.6 


^ 


>.78 
s(s+9ol) 


B. 


A  = 


1 

0 


s+31.5 
0 


0 
0 

1 

850 


>-78  ^ 
s(s+9ol) 


C. 


s+17.6       -'■ 

s^  +   58.2s  +  lOOls^  +  50^0s^  +  60O9  X  lO^Ki s  + 
p 


6.09  X  lO^Ko 


s^(s+9ol)(s+17o6)(s+31.5) 

Char.Eqn.  =  s^+58.2s  +1001s^+50^0s^+6.09  x  10^Kj_s+6o09  x  10^2=0 
(Order  of  CD  a  )  =  (3,11) 

Phase  II;   Relating  Specifications  to  Coefficients  of  Model 
Characteristic  Equations. 

A.   Statis  Specifications g 

1.  Steady  State  Position  Error  in  response  to  a 
step  =  0 
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a^  =  Limit  ^G(iDAll^.CD^  =  6.09  x  10% 

As  before,  the  value  of  a^  is  independent 

of  this  specification. 

Following  Error  in  response  to  a  ramp  =  0 

Error 


Limit  1 
s  ->  0  s 


1  -  ^i  A- 


=  Limit   s^  +   58. 2s^  +  lOOls^  +  50^0s 
s  ->•  0         Char.  Eqn. 

Therefore,  the  values  of  the  coefficients  are 

independent  of  this  specification, 

3.   Response  to  Load  Torque  <  10~  rad/ft.lb. 


10 


10 


>  Limit 
s 


^    [Char. Eqn.   J 


>  Limit   (.855)(^.78)(s)(s-H7.6)(s+^1.5) 
s  -+  0         Char.  Eqn. 


B, 


10  ^  >  0 

Therefore,  the  specification  on  load  torque 

disturbance  is  met. 
Dynamic  Specifications: 
Same  C  and  co  are  required. 


C.   The  same  dominance  factors  are  acceptable. 
Phase  Ills  Determining  the  Effect  of  Possible  Compensators. 

The  (Order  of  G^)  for  the  G^  feedback  position  is  changed 

■fco  (1,-1)  by  the  addition  of  the  parallel  integration.   The 

(Order  of  CD/x  )  is  changed  to  (3,11), 

"(2,1)  (2,11)  (1,11)  (0,1) 

(0,1)  (2,11)  (1,11)  (0,1) 

(1,1)  (1,1)  (2,11)  (1,1) 

(2,1)  (2,1)  (1,1)  (2,1) 


New  Location 

Matrix 
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Phase  IV:   Choosing  the  Compensator. 

A.  Using  °i  =  3  and  a^  =  2,  the  minimum  fifth  order 
model  is: 

s^  +  35.53s^  +  i527s^  +  ^-ajioos^  +  M-6^,ooos  +  3,030,000=0 

Since  this  is  now  a  Type  II  system,  it  has  zero 
error  in  response  to  inputs  of  position  and  velocityo 
The  a-./a  ratio  represents  velocity  error  in  response 
to  an  acceleration  input. 

The  requirement  that  the  working  model  meet  static 
specifications  is  removed. 

B.  1.  Since  no  working  model  had  to  be  formed  to  con- 
sider the  static  specifications,  the  initial  comparison 
will  be  made  to  the  minimum  model.  With  any  easily 
realizable  values  for  K-j^  and  K2,  the  a^^,  a-j^  and  a^  of 
the  uncompensated  equation  are  above  those  of  the  minimum 
model.  Since  a^_  is  the  only  fixed  coefficient  above 

the  minimum  model,  it  will  be  used  as  a  satisfactory 
coefficient.   K^  and  K^   will  be  left  unspecified, 
2.  n-m-2=l  and  m=5-2-l=2  arbitrary  dominance 
factors.   Solving  the  parametyic  form  of  the  a)^ 
coefficient  of  the  oscillatory  working  model, 
with  ci-i=3  and  02"^ »  yields  § 

(2C-L  +  2a^C2  +  °-2^%   "  ^^'^ 

(1.2  +  6C2  +  2)11,1  =  58,2 
^2  =  »53^ 
This  non-dominant  damping  ratio  is  larger  than  needed, 
so  ^2  =  »3  can  be  used  to  resolve  for  a  larger  a^o 
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(1,2  +  1.8  +  ^^)\\.\   =  58.2 
a^  =  2o2^- 

Using  these  values  for  a's  and  C's,  the  remain- 
ing coefficients  can  be  found, 
a.  =  2820 


a-L  =  56^000 
a^  =  3^15000 

The  compensated  characteristic  equation  iss 

s^  +  58. 2s^  +  2820s^  +  5^600s^  +  56^000s  +  3^-15000  =  0 
C.   1.   The  coefficients  to  be  changed  are  a^,  o.^t   ^i 
and  a  .  a-,  and  a  could  be  changed  independently  of 
additional  compensation  by  adjusting  K-j^  and  K2,  but 
K2/K,  has  been  set  at  2  5.  Therefore  it  is  necessary 
that  compensation  provide  for  varying  a-,  and  a^  indepen- 
dently of  the  other.   In  order  to  keep  the  coefficients 
to  be  changed  successive,  Kp  will  be  determined.   From 
the  characteristic  equation  developed  in  Phase  I-C  and 
the  compensated  equation  developed  in  Phase  IV-C-2§ 

60O9  X  10^  K2  =  3Al5jOOO 

K2  =  .56 
Therefore,  K^  =  K2/25  =  .022^- 

The  "semi-compensated"  characteristic  equation  is  now? 

s^  +  58o2s^  +  lOOls^  +  50'+0s^  +  137,000s  +  3Al59000=0 
The  compensated  characteristic  equation  iss 

s^  +  58. 2s^  +  2820s^  +  5^,600s^  +  56^, 000s  +  3,^159000=0 
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Comparison  of  these  two  equations  gives  the  desired 
Coefficient  Indicator  to  be  (2,1). 

2.   The  compensator  should  be  the  simplest  one 
which  has  sufficient  adjustable  parameters  to  change 
the  required  number  of  coefficients.   The  compensator 
must  have  three  adjustable  parameters  since  there  are 
three  coefficients  to  be  raised. 

The  compensator  must  be  in  a  feedback  path  since 
feed  forward  compensators  reduce  the  size  of  the  coef- 
ficients they  effect.   Comparison  of  the  desired  Coef- 
ficient Indicator  and  the  Location  Matrix  indicates 
that  feedback  position  G^^,  requires  the  simplest  com- 
pensator with  necessary  parameters,  one  with  order 
number  (2,0).   This  can  be  achieved  with  a  transfer 
function  of  the  form  K(s+a)(s+b)  which  will  arise 
from  the  following! 


to  b^ 


K-^s' 


K^_s 


from  G 


dD 


PF  =  K3(s^+Vk  s+^5/k  ) 

3«  A,  =  A-GiAj 

Z\ji  =  s^+58.Zs^+1001s^+50^-0s^+1^7,000s-H^Al5,000 
s  (s+9«l)(s+17c6)(s+33.o5) 


.  (^o781(l500)<850)K^(s^-HVK:>3-f^5/K^) 
s(s+9.^1)Cs+17o6)(s+31o5) 
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la  C1D--1M-  V  s(s+9ol)/l  s  s+17c6   s+31o5J 


G^  G 


la^dlAl^  ^   1^7»000(s+25) 

^m       s^+58o2s^+(6o09xlO^K^+1001)s^— .^ 


^  +(60 


09x10^X^^+50^0)  s^+ (6  o09xlo\^)s+3Al5, 000 


Equating  coefficients  of  the  two  forms  of  the  compen- 
sated equations  give  the  following  relationships; 

a^  =  60O9  X  10^  K^  +  1001  =  2820   or  K^  =  ,000299 

a2  =  6c09  X  10^  K)^  +  50^0  =  5^,600  or  K^_   =  0OO813 


a-^  =  6,09  X  10  K^  =  56^,000 


or  K, 


0O925 


8.^  Checking  the  Solution 

The  compensated  system  has  the  following  block  diagrams 


Loati  Torque 


s 


,022^ 


0855 


^ 


1500 

s+3io5 


-o- 


850 


s+17o6 


^■ 


^o78 
s(s+9ol 


f-^' 


,008135 


9r 


000299s' 


0925 


A. 


(o022^-+.56)    1 


1500 
s+31o5 


0  - 
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>o78  ^, 
s(s+9ol) 

[  ^o78 

(o299xlO"^s^ 

s(s+9ol)J 

+0813x10"^ S+0O925) 

0 

850 


s+17o6 


A   =  s^-^58c2s^-^2820s^-^5^,600s^-^56^,000s  +   ^^^15.000 
™  s^(s+9ol)(s+31.5)(s+17o6) 


For  an  input  of  a  unit  step  to  node  a  and  output  from 

block  G,j^,  the  system  performance  function  is§ 

1^7.000(s-H25) ^ 

^^  ^  s(s+2^.9)(s^+20s+1109)(s^+13o3s+123o2) 

The  response  of  the  system  in  time  domain  can  be  found 
by  the  method  of  evaluating  residues  to  bes 

©^^^(t)  =  l-^o0011e'^^°^^+ol37e"'^°^cos(31o8t-30°o6) 

+Ioi+8e"^°^^^cos(8o88t»139°o5) 

It  is  apparent  that  the  contribution  of  the  real  root  is 
negligible,  being  only  on  the  order  of  1  /  at  its  greatest o 
The  contribution  of  the  non-dominant  oscillatory  roots,  though 
somewhat  larger,  is  still  less  than  one-tenth  of  the  contri"= 
bution  of  the  dominant  roots o   It  decays  quickly  and  by  t  =  o2 
adds  only  about  2°/q  to  the  value  of  the  output o  At  t  =  06  sec, 
the  value  of  the  output  is  about  o98,  a  2°/  deviation  from 
final  value o   A  1°/^  deviation  was  specif iedo   At  t  =  o27, 


81 


peak  overshoot  occurs  and  equals  lo2^5  compared  to  the  speci- 
fied 1.20.   Failure  to  meet  specifications  exactly  in  these 
two  instances  is  attributable  to  placing  the  zero  too  close 
to  the  origin.   This  indicates  a  need  for  further  study  of 
the  effect  of  zeros  upon  the  transient  response,  or  for  a 
more  conservative  choice  of  a's  in  connection  with  zeros o 

The  solution  obviously  meets  static  specifications 
because  of  the  introduction  of  the  feed  forward  integration 
compensation. 
8.5  Summary  of  the  Application  of  the  Technique  to  this  Problem 

The  solution  of  the  problem  proceeded  along  straight- 
forward lines  until  it  was  discovered  that  a  Position  Control 
System  could  not  provide  acceptable  following  error.   It  was 
then  necessary  to  stipulate  compensation  which  would  give  the 
system  Velocity  Control  characteristics.  With  this  compen- 
sation inserted,  the  solution  was  then  recommenced.   The 
necessity  of  another  compensator  became  apparent  and  its 
characteristics  were  determined. 

It  is  interesting  to  note  the  similarity  between  the 
compensation  determined  by  this  solution  and  that  provided 
by  Chestnut  and  Mayer.    Both  involve  feedback  of  first 

and  second  derivatives  of  the  output o   However,  this  solution 

2 
replaces  their  elaborate  filter,   K5  ''s+a)        ,  with 

(s+b)(s+c)(s+d) 

feed  forward  through  the  function  K  and  feedback  of  a  signal 

s 

proportional  to  the  output  quantity. 
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CHAPTER  9 

CONCLUSIONS  AND  RECOMMENDATIONS 

9.1  Conclusions o  Application  of  step-by-step  procedures 
developed  in  this  thesis  resulted  in  successfully  designing 
compensation  for  a  system  which,  though  not  overly  complex, 
was  subject  to  somewhat  severe  specif icationso  Alternate 
methods  of  compensation  design  may  have  resulted  in  equally 
simple  comepnsation  but,  by  their  trial-and -error  nature, 
must  be  assumed  to  be  less  straightforward »   It  it  there- 
fore concluded  that  the  proposed  compensation  technique  is 
valid  and  presents  the  designer  with  a  powerful  and  labor- 
saving  tool. 

9 .2  Rpcnmmpndations »  Several  areas  for  further  study  present 
themselves  upon  consideration  of  the  results  of  this  thesis 
and  of  the  steps  by  which  these  results  were  obtained » 
Specific  areas  requiring  further  work  are: 

ao  More  exact  determination  of  the  deviation  from 
second-order  response  caused  by  the  existence  of 
zeros  and  additional  roots « 

bo   Careful  consideration  of  the  step-by-step  procedure 
to  eliminate  possible  duplication  of  effort  when  order 
raising  compensation  is  usedo 
Areas  of  research  which  should  provide  beneficial  results  ares 
a.   Examination  of  the  means  by  which  specifications  can 
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be  related  to  dominant  third-order  response. 

b.  Consideration  of  creating  non-dominant  oscillatory 
modes  at  harmonics  of  the  dominant  mode,  in  order  to 
better  reproduce  the  wave  form  of  the  input. 

c.  Examination  of  the  means  by  which  this  method  may 
be  extended  to  include  treatment  of  other  than  posi- 
tional servomechanisms. 


Qh 


APPENDIX  A 

A  DETAILED  EXAMINATION  OF  THE  EFFECT  OF  a  ON 

THE  DEVIATION  OF  A  THIRD  ORDER  RESPONSE  FROM  THAT 

OF  A  SECOND  ORDER 

To  examine  the  effect  of  a,  on  the  deviation  of  the  third 
order  response  from  the  second  order,  consider  the  variation 
of  the  two  systems  if  a  unit  step  is  applied  to  each. 

Second  Order:   Out(s)  =  2 ~ 2~  (A.l) 

s(s  +2Ca3„s+a)^    ) 
n       n 

Out(t)   =  -^  + \ 2  e"^^n^   sinCvT^^co  t+lf.) 

n  n 

(A, 2) 

where  ^^=  Phase   lag   of   ideal    (2nd  order)   response  = 
2 


-  arctan         r 

Third  Order:      Out(s)   =  ^'-^ 2~  (^'3) 

s(s+aTCO„)(s  +2Ca3„s+w     ) 
in  n       n 

where   ^j^^  is   included  to  make   steady  state  values  equal 

•acc_t 


Out(t) 


e n 


co^"       co^^[a^-2?:d+l] 


-Ceo  t 


1 n     ^c 


<«„^[(l-e){a2-2Ca*l)]V2  (AA) 
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where  lt^=  Phase  lag  of  actual  (3rd  order)  response  = 
a 

-  arctan  ^/l-C   -  arctan  n/1-C   . 
-C  a< 

Define  deviation  as  follows: 

D  =  actual  -  Ideal  (A. 5) 

ideal 

Substitution  of^'equations  (A. 2)  and  (AA)  into  equation 
(A. 5)  yields: 


p2  ^-iaO^t   ^    ^    3,^(^^2  ^^,^^^^   sinCv^^.t.*,: 
2.2Ca.l]  [a2-2U.;] 


T<^  e'^'^n'   +  sin(v^]7C^  co^t  +  f^) 


(A. 6) 


This  deviation  is  a  function  of  t,a,  and  C  of  the  dominant 
roots.   Ideally  we  would  like  to. find  the  maximum  deviation, 
no  matter  when  it  occurs.   This,  however,  is  a  great  under- 
taking which  is  beyond  the  scope  of  this  paper.   To  simplify 
the  problem  we  have  picked  the  time  corresponding  to  the  peak 


overs] 

loot  of 

the    i 

second  order  response: 

\- 

.1^ 

^n 

•    ■"^-  •  ^ 

Substituting 

this 

time 

into  equation    (A. 

.6): 

^ 

/nc^e 

-(a-c; 

'Vp 

a   sin(~»-arctan  -^"^ 

D  =L 

a^-2Ca+ll 

v^ 

<'  e- 

^'^n'^V  +  1 

(A, 7) 
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.v/1^, «-l 


but  Sin  (f+arctarr-^^^)  =  — ^  ^'^ t-tx 

2        ^-^      [a^-2Ca+l)-^/2 


(A. 8) 


This  equation  is  plotted  in  figure  A.l.   Figure  A. 2  represents 
the  result  of  similar  consideration  of  the  effect  of  a  zero 
at  ao)  on  a  second  order  response. 

Examination  of  these  curves  will  yield  an  appropriate 
value  of  a  to  insure  the  desired  degree  of  dominance  of  the 
second  order  response. 
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APPENDIX  B 
THE  MODEL  CHARACTERISTIC  EQUATIONS 


B-1 

Second-Order  Equations 

s^  -H  a^s  +  a^  =  0 

Minimum  Model: 

^1  = 

2C«„ 

2 

Working  Model: 

^1  = 

2Ca>„ 

2 

B-2 

Third-Order  Eauations 

^      2 

r  +  a^s  +  a, s  +  a^  =  0 
2      1     o 


Minimum  Model: 


a2  =  (2C+a-L)co^ 
a^  =  (l+2Ca^)a3^^ 


a^  =  (a,  )co. 


3 


o  In 
Working  Model: 

a^  =  (l+2Ca^)co^^ 

a^  =  (a  )a)^-^ 

o  In 


9b 


B-3  Fourth-Order  Equations 

h  ^  2 

s  +  a^s-'  +  a^s  +  an  s  +  a^  =  0 
3      2      1     o 


Minimum 

Model 

«3 

= 

2Cco„ 

«2 

= 

(l^a^2)<o/ 

^1 

= 

(2Ca^2)^^3 

^o 

= 

Working  Model  (Non-dominant  roots  complex) 
^3  =  (2Ci+2aiC2)a3„ 

«2  =  (l"^"l'^lV°l')%^ 
^1  =  (2a^C2+2t:ia^2)^^3 

^o  =  ^°1  '"n 

Working  Model  (Non-dominant  roots  real) 
a^  =  (2C+P-j^)a)j^ 

a2  =  (l+2CPi+P2^'^n^ 
^1  =  (Pl^^CP^)^)^^ 


%  =  (P2)-n 


where 


"r"2 


2  =  °l"2 


B-^  Fifth-Order  Equations 

s^  +  a^^s  +  a^s  +  ^2^     "^  ^1^  "'" 
Minimum  Model: 

\   =  (2Ci+a2)a>^ 
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B-^  Fifth  Order  Equations  (cont'd) 
Minimum  Model  (cont'd) 


^2  =  (2Cia^2^a2-Ha^2^2)a)^3 

Working  Model  (One  pair  non-dominant  roots  complex) 

a^  =  (2C;L+2a^C2+a2^"'n 

a^  =  (2aiC2+2Cia^^+a2^^a^a2CiC2+a^2  ^2^0)^3 

2  2  h 

a^  =  (a^  +20^0^^2+20^  °2^1^"% 

%  =  ^^1  "2K 
Working  Model  (Non-dominant  roots  real) 

a^  =  (l+2Z-^»^+»^)m^^ 


^2  =  (^*2!:jP2+P2)<o„3 
^1  =  ^h^^'^lh^'^n'' 


%  =    (^3^% 


5 


where 


°1  *  "2  *  °3 


'2  =  '"l'"2  "■   °2"3  *  "S"! 


^l'"2°3 
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B-5  Sixth-Order  Equations; 

s  +  a^s''  +  a^^s  +  a-.s"^  +  a2S  +  a-j^s  +  a^  =  0 
Minimum  Model: 

2   2   2 

^2  =  ^°1  -"^2  ""^1  °2  ^^n 

^1  =  ^^h\^2^^^n 

,   2  2,   6 
%  =  ^^1  "2  ^^n 

Working  Model  (Both  pairs  non-dominant  roots  complex) 

a^  =  (2C-L+2a^C2+2a2C3)a)^ 

a^  =  <l-^^^1^1^2'"^l^'''^^1^3V^^2S^lV^2^^^n^ 
^3  "  (.2a^r:^+2r:^a^^+2a^r:^-^8a^a^r:^r:2^^+2r:^a^^a^ 


^2r:^-2^^2a^a^\)c^^^ 


^2   "  (<ii^+^ai°^2^2^3"^^''l^'^2^1^3''^2''"^^^1^2''l''2' 


f      2     2.      6 
%   =    ^^1    "2  )^n 

Working  Model  (One  pair  non-dominant  roots  real) 


a^  =  (2C-L+2a^C2+Pi)a)j 


1^2^1^^n 

i^.  =  (l-Hl+a^CiC2+a^2+2C-Lp-L+2C2^iPi+P2^^n^ 
^3  "  (20^^2+2^3^ a^2+p^+Ifa^p^C^C2+a-^^P^+2C-LP2 

^2-1^2^2)^n^ 
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B-5  sixth-Order  Equations  (cont'd); 

Working  Model  (One  pair  non-dominant  roots  real) 

/■  2.  ,   6 


where 


1  =  V°3 


2  =  "a-'a 


Working  Model  (All  non -dominant  roots  real) 
aj  =  {2Ci-Pi)<c„ 
a^  =  (l+2':iPi^P2)a>„2 

«2  =  (^2*2CiP3+^)<»„'' 
(Pi.) 


where 


^o  -  ^^i+^"^n 


1  =  °l-'°2-'°3"'% 


^2   =   "lV°lV°lV°2  V2V"3\ 
^3  "   U^2/^3"'°1^2V°2°3V°1^3^ 
Pi,  =  °1^2"3^ 
B-6     Seventh-Order  Equations; 

s     +  a^s     +  a^s^  +  ai^s     +  a^s^  +  a^s^  +  a-^^s  +  a^  =   0 


Minimum  Models 


9^ 


Working   Model    (Some 
\  = 


^2   = 


^1  = 


2Ci+a3)co^ 

°l'^''2'*"l'°2^^2Cia^2„^^2C,«/a3)a>„'* 

2':i"l''-2'^-l'-3^"2'V"l^°2'°3^"n^ 

2      2^.,        2      2,6 
°1    '"2  -"^Ci"!    -a    "3^% 

2      2,7 
"l    "2    "iK 

roots  complex) 

V°^^n 

V°^5^^n 

V"^2^^n 


■^^■'"^At  )a)„ 
o       In 


o     n 


where  A  refers  to  coefficients  of  same  type  model  in 

sixth  order  eqns« 

Working  Model  (All  non-dominant  roots  real) 


^3  =  (f2^2!:iP3+Pi,)(o„ 
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where 


^2  =  (P3^2'^lN*^5>»n 


1  =  "l^  WV°5 

'3  =  °l''2V°l''2V°i'"2'"5*°2°3V"'2''3°5*°3^''5 


*"'l'"3V°i'"3''5*°l\''5'"'2%°5 
1,  =  ■'i°2''3V'"i'"2°3''5*''2''3\''5*°l''3\"'5'"°l°2\''5 
'5  =  °l°2°3V5 
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